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We take the first steps towards identifying the hydrodynamics of group field theories (GFTs) and 
relating this hydrodynamic regime to classical geometrodynamics of continuum space. We apply to 
GFT mean field theory techniques borrowed from the theory of Bose condensates, alongside standard 
GFT and spin foam techniques. The mean field configuration we study is, in turn, obtained from 
loop quantum gravity coherent states. We work in the context of 2d and 3d GFT models, in 
euclidean signature, both ordinary and colored, as examples of a procedure that has a more general 
validity. We also extract the effective dynamics of the system around the mean field configurations, 
^— ^ ' and discuss the role of GFT symmetries in going from microscopic to effective dynamics. In the 

process, we obtain additional insights on the GFT formalism itself. 

Introduction 

00 

Many approaches to Quantum Gravity have been developed and flourished in recent years [l| . Some of them build 
' ^j'" directly upon traditional quantization strategies (Dirac canonical programme, covariant perturbative quantization, 
sum-over histories or path integral quantization) applied to the gravitational field; examples of this are canonical loop 
quantum gravity and spin foam models in the continuum formulation 0, Q or the Asymptotic safety programme. 
(3J[)! Others rests on traditional quantization methods, usually path integrals, but introduce new ingredients or starting 
"~ assumptions, for example discrete space(time) structures interpreted either as mere regularization tools or as somehow 
| "physical" ; examples are simplicial approaches like quantum Regge calculus N|| and (causal) dynamical triangulations 
[5( , or spin foam models based on simplicial lattices [fj|, [Zj ■ String theory [8( is another example, as it starts with a 
^\ ' conventional covariant perturbative quantization strategy as applied to strings, and thus to their graviton excitations, 
but introduces a pletora of additional structures along the way, including extra dimensions, supersymmetry, and of 
course the abandonment of the point-like and thus strictly local nature of physical systems and their interactions. 
Other approaches still take their very start from some new radical hypothesis about the microscopic structure of 
f"^ quantum space and its constituents, variously inspired or motivated from General Relativity and geometrodynamics 
7—1 ! 1 or from more philosophical or mathematical considerations [1(1 HH ■ 

This paper deals with the approach to Quantum Gravity usually called Group Field Theory [l2l - fl5l ]. It is best 
understood as a generalization of matrix models for 2d quantum gravity [l6| . Matrix models, defined in terms 
, of pre-geometric, "space-free" constituents (matrices) with a dynamics dictated by purely combinatorial and discrete 
(pre-)geometric considerations, generate a sum over 2-dimensional simplicial complexes in their perturbative Feynman 
expansion, later to become a sum over smooth random surfaces in an appropriate continuum limit (understood in 
terms of a phase transition of the pre-geometric system) , defining a path integral for 2d quantum gravity. In a similar 
way group field theory models are defined in terms of "pre-geometric", "space- free" fundamental entities, fields over 
several copies of group manifolds or their Lie algebras, with a dynamics again dictated by purely combinatorial, 
discrete (pre-) geometric inputs, and generate d-dimensonal simplicial complexes in their Feynman expansion. Their 
historic origin from work in simplicial quantum gravity, state sum models and, more recently, in loop quantum gravity 
and spin foam models is manifest in the very mathematical structures they are based on. More precisely, while 
the combinatorics of field arguments in the GFT action is dictated, as in matrix models, by the requirement that 
the corresponding Feynman diagrams have the combinatorial structure of simplicial complexes of the appropriate 
dimension, the arguments of the field themselves and the specific form of the GFT action are motivated by results 
in loop quantum gravity and discrete gravity (and non-commutative geometry). As a result the GFT Feynman 
amplitudes take generically the form of simplicial gravity path integrals [l7| or, equivalently, spin foam models fl8j |. 
i.e. sum-over-histories of spin networks. The GFT field itself can be interpreted as a second quantized simplex or 
spin network vertex wavefunction. We refer to the reviews [l2l - fll| for more details. 
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Most of the above approaches, including some of those starting as straightforward quantizations of continuum 
General Relativity, end up identifying purely algebraic and discrete structures as fundamental building blocks (states) 
of quantum space, each characterized by a finite number of quantum degrees of freedom. However convincing the 
proposed microscopic dynamics for such fundamental building blocks is, the main task that all such approaches have 
to fulfill becomes that of showing how the algebraic, discrete quantum degrees of freedom give rise to a continuum 
geometric description of spacetime with its dynamics governed by General Relativity, in some approximation. 

Various solution to this problem have been proposed. These range from coherent states in canonical loop gravity 
[l9j to statistical field theory methods in simplicial quantum gravity H, Q and in other discrete pre-geometric models 
of space [20| | . to ideas coming from quantum information theory [211 ] . Unless a suitable interpretation in terms of 
continuum geometries is found for both the discrete quantum structures and their microscopic dynamics as such (for 
attempts in this direction see H24i3])j a continuum space is likely to be understood as built out of a very large 
number of microscopic quantum building blocks, and the continuum limit of the microscopic dynamics arises in a 
thermodynamic limit, and is thus possibly subject to critical phenomena. This is the idea at the root of simplicial 
gravity approaches, hence the use of statistical field theory methods, but also the origin of the very difficulties involved 
in tackling the problem of the continuum, since one has necessarily to deal with very large numbers of quantum degrees 
of freedom and their collective behaviour. The same difficulty is encountered in the context of loop quantum gravity 
coherent states and spin foam models, even more so if dynamical spatial graphs or spacetime lattices are involved (as 
the graph-changing nature of the Hamiltonian constraint operator in loop quantum gravity would seem to require) . 

At this point, group field theories enter the stage as the framework in which, on the one hand, the very same 
fundamental structures identified by other approaches as building blocks of space (simplices or spin networks) can 
be treated in a unified way, with all the techniques already developed in those approaches. On the other hand, new 
techniques become available, thanks to the GFT embedding. In particular, GFT become the natural formalism for 
studying the physics of a large number of the same quanta of space identified in other approaches, because it is a 
second quantized formalism of the same, well-adopted to the study of many-particle systems. Let us clarify this point 
further, as it is central to the work presented in this paper. 

Once understood that the same building blocks of space used in other approaches arise as quanta of the GFT field, 
and that the GFT action defines a microscopic dynamics for them, then the study of their effective dynamics (classical 
or quantum) in a thermodynamic limit becomes a problem in statistical GFT. Quantum space is then interpreted, 
at least at the level of analogy, as a sort of exotic condensed matter system, a condensate/fluid made out of large 
numbers of GFT quanta, to be studied with tools and ideas from condensed matter theory. The continuum and 
semi-classical approximation that leads from microscopic quantum structures to smooth continuum space(time) can 
then be seen as analogous to the hydrodynamic approximation that leads from the microscopic quantum description of 
a few atoms system to the hydrodynamic description of a classical fluid made out of large numbers of the same atoms. 
The GFT itself becomes then the quantum gravity analog of the microscopic quantum field theory of non-relativistic 
atoms that underlies any condensed matter system, and any (quantum or classical) fluid. Classical geometrodynamics 
(including General Relativity) should then arise, in a way to be understood, from GFT hydrodynamics, in some 
appropriate phase of the theory, also to be identified. It is unlikely, in fact, that the infinite number of GFT degrees 
of freedom that are necessary to obtain a smooth spacetime and its relativistic dynamics will organize themselves in 
a single macroscopic phase, and thus that no understanding of GFT phase transitions will be needed to understand 
the emergence of continuum classical spacetime from GFT. This perspective has been argued for in [25| . 

The idea of spacetime as a fluid/condensate has of course been put forward repeatedly (26|, and is somehow 
the conceptual underpinning of condensed matter analog gravity models [13, HH , together with the idea of General 
Relativity as the emergent hydrodynamics or thermodynamics of microscopic pre-geometric building blocks (291433 ] . 
The body of this work (in particular the one on quantum fluids 1 ), we believe, represents an important guide for the 
extraction of geometrodynamics from GFT. Conversely, we also believe that GFTs can represent a candidate for the 
microscopic description of the "atoms of quantum space" that all this work somehow hints at. 

Once accepted, at least provisionally, the above perspective, the task is to identify the best strategy and mathemat- 
ical tools to tackle the problem of the continuum in a GFT context. Being a problem in statistical GFT, obvio usly 
the renormalization group is a key asset. Indeed a programme of GFT renormalization has recently started [371 ] . 
which has extracting the continuum limit of GFT as its main physical goal, beyond the many mathematical insights 
it is providing ]47] . Another key method for extracting information about the phase structure of a condensed matter 
system, and for the extraction of the corresponding effective dynamics, is mean field theory. In particular, mean field 
theory is the most direct route from the microscopic atomic dynamics to the effective hydrodynamics in the case of 
quantum fluids and Bose condensates, which are the systems we tentatively use as paradigms of what the emergence 



1 For concrete implementations of these ideas within idealized BEC models, and extended discussions, see [35, 36]. 
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of spacetime from quantum gravity may entail. 

The application of mean field theory requires first of all the identification of a candidate macroscopic configuration of 
the system (and associated collective variables), whose dynamics encodes the collective dynamics of the many- particles 
constituting it. Also, the fluctuations around such new vacuum become the relevant (once more, collective) degrees 
of freedom at the macroscopic level, with an associated dynamics in general very different from the microscopic one 
governing the underlying atoms. All this works under the assumption that this new vacuum will be very different 
from the microscopic quantum vacuum of the system (the Fock vacuum), and that the physics one is looking for at 
the macroscopic level is best understood close to the assumed macroscopic configuration, rather than the microscopic 
Fock vacuum. This a particularly reasonable assumption in the GFT context, where the vacuum around which the 
usual perturbative expansion takes place is interpreted physically as a "no space" state [l2l - [l5j . in which no space 
structure at all exists, neither topological nor geometrical. The same reasoning suggests that the usual spin foam 
models or the equivalent simplicial path integrals are not the most convenient definitions of the dynamics of the degrees 
of freedom that correspond to a smooth s pac etime and its geometry, because they arise as Feynman amplitudes of 
GFT models around the no-space vacuum |25j |. It suggests to look instead for an effective dynamics of perturbations 
around a different vacuum. The identification of the relevant macroscopic vacuum and of the corresponding dynamics 
is however a highly non-trivial task, in condensed matter and, even more, in our quantum gravity context. The best 
one can do, usually, is to proceed by educated guesses, and then test the resulting hydrodynamic theory against 
observation. In quantum gravity we cannot (yet) make use of experimental inputs, but we can still proceed using 
intuition and the large amount of theoretical ideas accumulated so far, and study the formal aspects of the effective 
theory obtained as result of our hypothesis. 

This is what we do in this, largely exploratory, paper. We move the first steps towards establishing a hydrodynamic 
limit of GFTs and in relating this hydrodynamic regime to classical geometrodynamics of continuum space. First 
of all, this requires the identification of a candidate (non-perturbative) macroscopic vacuum. Here we start from 
the results obtained in the context of semiclassical loop quantum gravity, in particular on LQG coherent states. We 
identify, from a simple analysis of LQG coherent states, a candidate coherent state wave function associated to vertices 
of spin networks. This coherent state wave function will depend on parameters that can be interpreted in geometric 
terms, following again the LQG results, and that play the role of order parameters in the GFT context. We identify 
this candidate vacuum state in section III, together with its geometric interpretation. Then, we re-interpret the same 
wave function as a classical GFT field, and use it as candidate mean field configuration for the GFT dynamics. We 
do so in section IV and V, considering in detail the simpler cases of GFT models for quantum BF theory in 2d and 
for 3d quantum gravity, in euclidean signature, and for both ordinary and colored GFTs. We extract the mean field 
theory equations for the order parameters, as resulting from the GFT dynamics. These hydrodynamic equations 
have a geometrodynamic interpretation, that we elucidate to some extent. We are not able, at this stage, to give a 
complete and transparent geometric re- writing of these equations, nor to obtain a precise mapping of these equations 
with those of General Relativity. This is perhaps the main limitation of our results. Still, we believe we open up 
a new interesting direction for further investigations. Having obtained the mean field or hydrodynamic equations 
for the (geometric) order parameters characterizing this GFT vacuum, we move on to extract the effective dynamics 
for perturbations around it, in section VI. The main purpose, here, is first of all to illustrate the general procedure 
and its main features; second, it is meant to show how, concretely, the effective dynamics could, on the one hand, 
differ from the microscopic one and, on the other hand, depend on it, in the spin foam formulation. A more detailed 
analysis that the one we perform, once more, is certainly needed, but most important, what is needed is a careful 
study of the physical interpretation of the perturbation field itself, and of the degrees of freedom it carries. At the 
present stage, it is unclear if these should be interpreted as matter fields 2 living on the geometric background defined 
by the mean GFT field, or if they have a geometric interpretation as well, so that, for example, the full geometric 
character of the effective GFT hydrodynamics is to be looked for in the coupled equations for order parameters and 
quasi-particle fluctuations around the mean field. We leave this issue for future work. Last, in section VII, we discuss 
the issue of GFT symmetries, in particular diffeomorphisms, from the point of view of the application of mean field 
theory methods. A summary of our results, and some outlook, are provided in the last section. 

To clarify the path we follow and the conceptual framework for the results we present from section III onwards, it 
will be useful to briefly review how the mean field theory method is applied in a very well understood physical system, 
namely Bose-Einstein condensates, to extract the effective hydrodynamics of the condensate. We do so in section II, 
where we also point out the main differences (as well as the similarities) between our approach in the GFT context 



2 Recent work in the GFT context l38t |39H has applied a similar strategy to the one we follow in this paper, and can also be understood 
from a condensed matter analogy [401 ] - The purpose of such work was, in fact, the extraction of effective (non-commutative) matter 
field theories as perturbations around exact solutions of the GFT equations of motion. No effective dynamics for the background GFT 
configurations chosen was investigated. 



4 



and the case of Bose condensates, together with its limitations. 



I. MEAN FIELD THEORY AND EFFECTIVE HYDRODYNAMICS IN BOSE CONDENSATES: 
ANALOGIES AND DIFFERENCES WITH THE GFT CASE 

The paradigm we will follow is mean field theory applied to the study of Bose-Einstein condensates. For an extensive 
discussion on the subject, we refer to [ill ]. Here, we recall only the main features relevant for our purposes. 

We will focus on the particular case of Bose-Einstein condensation in dilute gases. For these systems, the most 
convenient treatment is provided by the second quantized formalism, based on quantum field operators 

*(x) = 5^a i « i (x), (1) 

i 

where i is an index labelling an orthonormal basis of (single particle) wavefunctions Uj, with 

w 4 (x)u*(x)d 3 x = Sij , (2) 

v 

and di are (bosonic) annihilation operators, obeying 

[ai,aj]=0; [fit, fit] = (3) 

These operators are associated to the creation/annihilation of fundamental particles (atoms) of the system. The 
typical Hamiltonian describing the dynamics of a dilute gas of bosons (in the local approximation) has the form: 

H = J *+(x) (-^- - M + V(x) + ~*t(x)*(x)) *(x) d 3 x, (4) 

where m is the mass of the bosons, V(x) is an external potential used to confine the system in a given region of space 
(V), k is a (positive) constant encoding the strength of two-particle interactions and /i is a chemical potential, playing 
the role of a Lagrange multiplier used to fix the total number of particles. The above Hamiltonian corresponds to the 
case of a single atomic species, but the formalism can be immediately extended (2BEC, spinor BEC, etc). 

The main problem is then the correct identification of the ground state. If the system were non-interacting, the 
ground state, below the critical temperature, would be the state in which all bosons occupy the same single particle 
ground (lowest energy) state (Bode condensate). The case of interacting particles is much more involved. However, 
one can make a guess at the macroscopic properties of the system by assuming that, below a certain temperature (to 
be computed), the condensation does take place even in the interacting case, and even for repulsive atomic interactions 
(positive k). The simplest way to implement this idea is to assume that the (interacting) ground state (G.S.) is a 
state in the Fock space, \G.S.), such that the operator $(x) acts on it as multiplication by a function, 

*(x)|G.S.) « tp(x)\G.S.) (5) 

Clearly, this kind of state is crucially different from the Fock vacuum (F.V.) \F.V.), for which if?(x)\F.V.) = 0. 
A simple example of state possessing this property is a (second quantized) coherent state. Consider the state: 

\ Zi ) =e-l z l 2 / 2 ex P (z4)|F.K) (6) 
It is immediate to see that this state is an eigenstate of the field operator ^(x): 

*(x)|zi) = z Ul (x)\z z ) = *(x)\G.S.}. (7) 
The content of the state is characterized by the expectation value of the number operator of each mode: 

Tij(i) = (zilajajlzi) = 6 lJ \z l \ 2 (8) 

Therefore, as a working hypothesis, one assumes that the ground state of the many body system is a coherent 
state, with a macroscopic occupation number of a suitably chosen single particle state: \G.S.) — \zi). Concretely, this 
means that one is working in a regime in which the field operator ^z(x) can be effectively separated as 

*(x)« V(x)I + x(x), (9) 
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where ip(x), often called the condensate wavefunction, encodes the information about the ground state, while x 
encodes the effect of deviations from the mean field ip. This splitting must then be introduced into the equations of 
motion for the field operator In turn, this equation of motion will become an equation for the classical field ip, 
including the effects of perturbations x- The logic is clear: the mean field ip must be determined self-consistently 
from the mean field equations. At the lowest order, neglecting the fluctuations x, this equation assumes the form 

d ft 2 V 2 

m— ip = — - — v - + + (io) 

at 2m 

known as the Gross-Pitaevski hydrodynamic equation. Its solution provides the description of the macroscopic 
properties of the condensate described by the wavefunction ip. It can be re- written in terms of two real fluid equations 
(Euler and continuity) for a perfect fluid of density p and velocity v such that 

iP = ^pe- ie , v<xV6. (11) 

It is in this sense that the mean field approximation gives immediately the hydrodynamic description of the system. 

It is important to keep in mind that this method must be self-consistent: we have started assuming that the ground 
state is a condensed state, in which many particles condense in the same single-particle state. However, it is only after 
that we have found the solution of the GP equation, with the appropriate boundary conditions, and we have established 
that this solution corresponds to a configuration in which a large fraction of the N bosons originally present in the 
system is occupying the condensed state 3 that we can say that our system condenses. Also, even if the method proves to 
be consistent, this mean field treatment is just an approximation that takes into account only the field configuration on 
which the coherent state is peaked. The spread of the coherent state itself is neglected (this corresponds to neglecting 
terms containing, e.g. (x 2 ))- I n general, more refined tools are needed (42J. Nonetheless, the rough description in 
terms of the condensate wavefunction allows to understand a number of interesting features of the physics of the 
system. From a perspective d la Landau, the condensate wave function is an order parameter describing the phase 
transition. It plays two roles. First of all, it contains the basic information about the hydrodynamical properties of 
the condensed fraction (density, velocity profile etc.). Second, it determines the symmetries and the general properties 
of the effective dynamics of perturbations around the condensed state itself. In fact, if we consider the dynamics of 
phonons, the elementary excitations above the condensate, their various properties (internal symmetries, dispersion 
relations and spacetime symmetries, etc.) are essentially controlled by the wavefunction ip. 

It is useful to make a comparison with the GFT case. In the BEC phenomenon the ground state being selected 
has the peculiarity of peaking the second quantized field operator on a given classical field configuration. From the 
point of view of the particle content, it is a superposition of many particle states with different occupation numbers. 
Finally, the mean field represents the order parameter that describes the condensation, and encodes the macroscopic 
hydrodynamic variables that effectively describe the semiclassical state. 

In the case of GFTs, the second quantization formalism, with a definition of a Fock space and a corresponding 
many particle interpretation is still not understood. Therefore, a proper GFT analogue of the BEC mechanism is not 
available cither. Consequently, in what follows, we will be still working in a first quantized scenario, corresponding 
to the kinematical Hilbert space of Loop Quantum Gravity, whose states appear as boundary data in GFT transition 
amplitudes, and therefore the analogy with the BEC case must be taken with care. The candidate ground state will 
be obtained from LQG coherent states, and correspond to wavefunctions for spin network vertices (the GFT quanta, 
i.e. their candidate 'atoms of quantum space') characterized by (peaked on) SL(2,C) group elements, defining points 
in the classical phase space appropriate for the gauge theories considered. The situation is similar to a BEC mean 
field approximation in which one assumes that the configuration uq (x) is a wavefunction peaked on a certain point of 
phase space for a single particle. This is clearly a much stronger semi-classicality condition than only requiring the 
quantum field operator acquires a nonzero vacuum expectation value, as in the true BEC case. 

Still, within the limits of our analysis, the two mean field approximations in BEC and GFT (as performed in 
this paper) share the same philosophy. The group elements where the peaks are situated play the role of order 
parameters, as the condensate wavefunction does. Also, they will be dynamically determined by equations that will 
be derived from the GFT equations of motion in the same way in which the Gross-Pitaevski equation is derived 
from the non-linear Schrodinger equation for second quantized field operator. As it has been discussed in details in 
the literature, the SL(2, C) group elements determine the mean value of certain geometric quantum operators, and 
correspond to classical geometric fields, while the GFT dynamics is expected to encode the full quantum dynamics of 



3 Explicitly we must check that f d 3 x\ip(x)\ 2 RJ N. 



FIG. 1: Entire oriented link from a to b, with associated group variable h a i, and peak group element g a b of the heat kernel. 



the microscopic degrees of freedom of quantum space, and the equations that we derive for them represent a form of 
geometrodynamics, here derived from the microscopic dynamics of quantum space as encoded in a GFT, in a mean 
field (hydrodynamic) approximation. 

II. FROM LQG COHERENT STATES TO AN ANSATZ FOR FOR A GFT BACKGROUND 

CONFIGURATION 

We now introduce briefly kinematical loop quantum gravity (LQG) coherent states, and then show how to extract 
from them a candidate background configuration to be used in the mean field approximation of GFT models. More 
precisely, we show how LQG coherent states can be obtained by convolution of vertex wave functions, also character- 
ized, as the LQG coherent states, by SL(2,C) group elements, and ready to be used as mean GFT fields, as we do in 
the following. 

A. Coherent states on SU(2) and vertex wavefunctions 

(Complexifier) coherent states were introduced within Loop Quantum Gravity in order to investigate the semiclas- 
sical limit of the theory. They are wavefunctions associated to graphs, which have the properties to be peaked around 
certain values of the classical phase space of the theory, parametrized by holonomies (group elements of SU(2)) and 
fluxes (elements of the Lie algebra su(2)). The expectation values of (polynomial functions of) these operators coin- 
cides with the classical values at the corresponding phase space point. Also, they minimize uncertainties for a specific 
set of kinematical observables and within a given class of quantum states. Given this nice kinematical properties, 
a key question becomes then to what extent this class of states satisfies the quantum dynamics of the theory, and 
whether the same dynamics preserve their semi-classical properties. While we will not address this question within 
the LQG canonical framework, we will in fact study the issue of whether the associated vertex wave functions, our 
candidate GFT mean field, solves the GFT equations of motions. 

For a complete discussion of their properties and their relevance for the semiclassical approximation within LQG, 
we refer to [191 ] . Here we summarize only their definition. 

Let 7 be a graph, and let g denote a set of group elements (SU(2)) associated to its links (one per link); these 
group elements should be considered as parameters defining the state. The coherent states on 7, in the connection 
representation, are functions from S\J(2) E ^ to the complex numbers, given by: 

¥[7, |(d.9) |yWI V[7,5;i]({5(e(l))^r 1 (e(0))}), (12) 

where h denotes the set of the group variables (one per link) , the integrations over the group associated to the vertices 
V(j) of the graph 7 enforce gauge invariance (Gauss constraint), and where the gauge- variant wavefunction, 

i>h,g;t](h)= J] Pt(h e gl l ), (13) 

is a product of heat kernels (pt) on the group manifold SU(2) [43[ with spread (diffusion time) t. Each heat kernel is 
a function of the group element h peaked on the element g with spread t. 

The same states can be rewritten in terms of functions associated to vertices, as we now show. For each link (ab) 
as in figure [TJ introduce the intermediate point P and double the associated group elements (see figure [2]) as: 

Kb = h^ha, g ab = g b ~ 1 g a , (14) 
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thus artificially doubling the number of group elements associated to the same link (this is not a physical doubling 
of variables, since the two group elements only enter the expression as the products h ahl g ab ). We then use the 
convolution property: 

P2t(9) = J dhpuigh-^pt^h); h+t 2 = 2t. (15) 

For simplicity (but without loss of generality 4 ) we will consider ti = t% = t. 

Neglecting for the moment the group averaging procedure enforcing gauge invariance at each vertex, we obtain: 



* = P2t(Kb9 a b) = P2t( h b 1}l a9a 1 9b) = P2t{ha9 a ^ 9b\ 1 ), (16) 

where we have used the property pt(hg) = p t (gh). Then, we can split the single link by using the convolution property 
of the heat kernel: 

dqab Pt(h a ga 1 Qab)Pt( < labgbK 1 ) = J dq ab Ptihaga 1 Qab) q~ b ) , (17) 

where we have used the property pt(.g _1 ) = Pt{g)- Notice the identical functional dependence on the various group 
elements in each of the two functions associated to each (oriented) semi-link. 

With an obvious change of notation (ab — > e, a — > e(0) and b — > e(l)), the complete function associated to the 
graph reads: 

#y,£2t](£)= f(dq)\ E ^\ J] PtiK^g^q^PtiK^g^q- 1 ). (18) 

J e€£( 7 ) 

This expression can be reorganized as 

#y,£;2t](ft) = j{dq)\ E ^\ J] #„(/K,^9^;^,2^Xl a ;».;^, ro „^< m J (19) 

where h V)i is the i— th element associated to the i— th link ending in the vertex v, m v is the valence of the vertex v, 
and: 

$v(gi;-;g mv ) = Y[p t (gi) (20) 

This vertex function is a function going from G m " — > K, thus can be interpreted as a GFT field. Correspondingly, the 
whole wavefunction can be rewritten in terms of a suitable operator defined in a group field theory 5 . Also the gauge 
invariance property fits in the same scheme. 
Imposing gauge invariance, we get 

&[ 7 ,g;2t](h)= j {d~g)\ v ^\{dq)\ E ^ \{ ^{{K^g^q^}) (21) 

i>ev(7) 

implying, in turn, 

&[ 7 ,g;2t](h)= j{dq)\ E ^\ \{ Mig^Ue'M), (22) 

where we have defined 

Mihi}) = J dg v $(g v hi), (23) 



4 Furthermore, the parameter t should be determined by the quantum dynamics {e.g. the Hamiltonian constraint), and a natural assump- 
tion is that the dynamics will select a value of t which will not depend on the portion of the particular graph we are considering, at 
least in the case in which we will not deal with too inhomogeneous states. 

5 Also the case in which the valence of the graph is not fixed could be accommodated within the GFT formalism, by using several fields. 
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FIG. 2: Splitting the link a to b. The group elements h a , hb, g a , gb refer to the semi-link from the vertices a,b to P, with the 
appropriate orientation. 

This analysis shows that the group averaged coherent states associated to a graph might be written in terms of a 
group field theory, with a suitable number of fields according to the valence properties of the graph itself. Therefore, 
we will use these particular field configurations within GFT as a form of test configurations. However, before doing 
this, we need to go one more step further in the definition of the semiclassical LQG states. 

In considering the coherent states for quantum gravity, in fact, one has to analytically continue the g a elements 
from SU(2) to SL(2,C) in order to peak the states on confi gura tion with given parallel transports and given triads 
per link. This analytic continuation exists and is unique fl9l |44| . 

In the case of complexified heat kernels (see appendix), the heat kernel itself is peaked around a SL(2,C) group 
variable. From here on, group elements denoted with a capital letter will correspond to SL(2, C) elements. The 
correct analytic continuation of the element associated to the link gives G a b — G b ~ 1 G a - In fact, using the isomorphism 
between SL(2, C) and T*SU(2), the cotangent bundle of SU(2), this SL(2, C) element is identified with a point in the 
phase space of LQG (or SU(2) BF theory) at that spacetime point, or, in the discrete setting, associated to any given 
element of the discrete spacetime lattice. In particular, considering cLgctin the single link, G SU(2), while g a b is 
replaced by G ab G SL(2,C). 

In order to discuss the implications of the decomposition, we need to introduce some notation. We recall the polar 
decomposition of the SL(2,C) matrices (making the correspondence with T*SU(2) explicit): 

Gab = gab exp(£ a6 ), G a = 9a exp(E a ), G b = gb exp(E b ), (24) 

where gx are SU(2) elements while Ex — E x Ui are su(2) matrices, with ui denoting the Pauli matrices. 

In the language of LQG, these SL(2,C) elements have a direct geometrical interpretation in terms of expectation 
values of holonomies and fluxes. More precisely, the SU(2) part will correspond to the expectation value of the 
holonomy on the link, while the E % will be interpreted as the corresponding classical value for the triad/flux operator 
associated to an elementary surface dual to the link and intersecting it at a single point. 

Using the above polar decomposition, we have 

g a bexp(E ab ) = exp(-S 6 )g^ 1 g exp(£' a ) = g b ~ 1 g a cxp(-E b ) exp(£ a ), (25) 

where we have defined 

E b = {9b l 9a)- l E b {g^ l g a ). (26) 

Let us discuss briefly the geometrical interpretation of the above group elements, considering two oriented semi-links 
coming out of two vertices being glued to form an entire link. 

If we follow the geometrical interpretation according to which the Lie algebra elements E a , Eb are (smeared) triads 
as seen from the vertices a, b, respectively, it is clear that they must be compared after appropriate transport on the 
midpoint at which the gluing takes place. Therefore, define (see figure [3]): 

E c = g a E a g-\ E d = g b E b g^ (27) 

In absence of gauge invariance, we should expect that, in order to have a sensible geometrical interpretation, these 
two triads must coincide modulo a sign (due to opposite orientations of the surfaces of integration) , 

E c = -E d . (28) 

Taking the exponential of this equation we obtain: 



g a exp(E a )g a 1 = g b exp(-E b )g b \ (29) 
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which leads to: 

Gag-^gbG^ 1 . (30) 

This last equation can be seen as a "geometricity constraint" on the data assigned to the vertices: if these constraints 
are satisfied, the links obtained by gluing the corresponding legs will be characterized by a very simple geometrical 
content, easily visualized in terms of gluing of cells (simplices, hypercubes or other kind of basic "chunks of spacetime" ) 
along common boundaries. This implies that: 

E a = -g^gbEtg^ga = -E b (31) 

whence: 

g ab exp(E ab ) = g^g a exp(2£ a ), (32) 

where the RHS is already a polar decomposition of the SL(2,C) element. Therefore the link, after gluing, is peaked 
on the parallel transport obtained by correctly composing the elementary parallel transports, while the triads will be 
uniquely determined by the triad as seen from one vertex, as we would expect 6 . 

Let us summarize briefly where we stand so far. The semiclassical canonical wavefunctions associated to a given 
graph have been written in terms of products of functions associated to the vertices. In turn, these vertex functions 
will have a dependence on certain group variables associated to the vertex itself and to the each leg attached to it. 
The combinatorics is such that, when assembled, the total wavefunction depends only on certain combinations of 
these variables, matching the original dependence on group variables. 

We have also seen that the possibilities for the peaks of the vertex functions are not always compatible with 
geometrical requirements. Rather, we have to impose some constraints, which will select, among all the possible 
assignments of SL(2, C) elements, the ones compatible with a consistent gluing of elementary geometric cells along 
common boundaries. 



It is interesting to note that, before imposing gauge invariance, we could be tempted to use these geometricity constraints to put more 
conditions on the data being assigned to the vertex functions. In particular, defining: 

Qx = G x g- 1 , (33) 

the geometricity constraints become: 

Qa = Q^ 1 (34) 

whenever we want to glue vertex a with vertex b through a particular leg. If we want to keep the geometric interpretation independently 
from the particular gluing pattern, i.e. if we want to keep the possibility of gluing every leg with any other leg, wc obtain a rather strong 
constraint. It is enough to consider these constraints in cycles: 

Qa = Ql 1 ) 

Qb 1 = Qc \^Qa = Qa 1 (35) 
Qc = Qa 1 J 

that in turn implies: 

Q a = ±1 E a = 0, (36) 

with the parallel transport left free. These constraints, then, would suggest that the only way in which the gluings arc consistent with a 
naive geometrical intuition independently of the gluing patter is that all the triads are vanishing. In fact, this is not the only possibility. 
One obvious way out is to include a constraint on the pattern used for the gluings, in order to avoid the ciclicity that we have exploited 
to derive our conclusion. 
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Additionally, notice that whatever the Lie algebra elements are, even if the vertex wavefunctions are peaked on 
SU(2) elements which are different from the identity, but they arc the same link per link, the wavefunction in terms 
of the graph is peaked on the flat connection. Similarly, after complexification and consider heat kernels peaked 
on SL(2, C) elements, whenever they are equal, the wavefunction associated to the link is peaked on the identity of 
SL(2,C) itself, i.e. not only on the flat connection, but also on zero triads. 

B. Candidate GFT mean field configuration 

Having decomposed the coherent wave function associated to a generic graph in terms of vertex wave functions, we 
can now re-interpret this special vertex wave function as a classical GFT field, and use it in that context. 
We will consider a compact Lie group G, and the group field theory defined from: 

4> : G m -> C (37) 

and a classical action S[4>] or equivalently classical field equations: 

where w is a label that identifies the various terms O w in the action. We limit our analysis to the case in which 
G = SU(2), for simplicity. We focus on the particular field configuration: 

/m 
dh\\p t {h l9 ^h) (39) 
i=l 

where g% are group elements (in SU(2) or in SL(2, C)) on which the heat kernels would be peaked if the integration over 
the group variable associated to the vertex would be absent; with such gauge invariance restriction being imposed, 
the function will be peaked instead on gauge invariant equivalence classes of the m group elements associated to the 
m links of the vertex. We ask the following questions: is £ a classical solution for some GFT? Is it an approximate 
solution for one of the cases we are interested in? What consequences, i.e. what restrictions on the parameters 
characterizing £, follow from demanding that it is in fact an approximate solution? In the corresponding regime of 
parameters, what is the effective dynamics for perturbations around such mean field configuration? In the rest of the 
paper we try to address these questions. 

To be more specific, we will consider several possibilities for the mean field £, in dependence of the possible choices 
that we have in realizing the GFT theory in terms of field content (with colors or not) and gluings (specific gluing 
patterns or not). These choices are summarized as follows: 

1. single field GFT 

a. all the links are peaked on the same SL(2, C) element; 

b. each link peaks on a generic SL(2, C) clement; 

2. colored GFT 

a. for each color, all the links are peaked on the same SL(2,C) element; 

b. for each color, each link peaks on a generic SL(2,C) element. 

While the situation in case (la.) leads necessarily to graph wavefunctions that are peaked on the identity in SL(2, C) 
for each complete link (i.e., after gluing), the others do not. According to the choice of group elements (which will 
be ultimately decided by the dynamics) and of the gluings (kinematical requirements on the particular GFT action) , 
the outcome might be a graph wavefunction which peaks, on each link, on nontrivial SL(2,C) elements. 

Our (longer term) goal in fact is a bit more ambitious than this. We want to establish a pattern (closely related to 
the self consistent mean field method that we have discussed in the previous sections) which applies to the problem of 
finding the physical content of certain group field theories, in a suitable macroscopic (hydrodynamic) approximation, 
and to relate it to possible dynamics of effective geometries. We will try to understand what is the form of the 
dynamical equations determining the mean field values (i.e. the peaks of the wavefunctions) once we make the 
educated ansatz that we should look for solutions in the form of heat kernel coherent states. As in the case of many 
body systems, the equations of motion for the full field will provide an effective equation for the mean values/peaks. 
These peaks will play the role of order parameters: not only they will determine the properties of the geometry on 
which the state is peaked, but also they will determine the shape and the symmetries of the theory obtained by 
expanding the GFT around the heat kernels obtained in this way. 
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III. 2-DIMENSIONAL CASE 



We start from a simple case, a GFT generating 2D simplicial complexes in its perturbative expansion, and corre- 
sponding to a quantization of topological BF theory with SU(2) as gauge group. Despite its simplicity, it will represent 
our case study, on which we will try to model higher dimensional cases. 

Before attacking the GFT as such, it is worth recalling some basic facts of the classical theory behind it (see [45[ 
for more details), that could be useful to give meaning to the later manipulations and results. The general form of 
the BF action is 

S = [ Tr (B A F(ui)) (40) 

J M 

where M. represents the spacetime D-manifold, B is a (D — 2)-form valued in the su(2) Lie algebra, w is a 1-form 
connection, also valued in the Lie algebra of the gauge group, and F(lj) is the corresponding 2-form curvature; the trace 
is taken in the fundamental representation of the algebra. The (spatial components of the) B field and the (spatial 
components of the) connection ui are canonically conjugate variables. At the discrete level, with the connection 1-form 
replaced by group elements representing its holonomies, and B field replaces by Lie algebra elements, one has that 
the classical phase space variables on each discrete element of the lattice is given by the cotangent bundle of SU(2), 
expressed in terms of the same variables. The classical equations of motion are 

d u (B)=0 F(u)=0 (41) 

enforcing the compatibility between connection u) and B field, and flatness of the connection. It is then clear that the 
classical theory is trivial, as it contains only flat configurations. At the quantum level, the B field acts as a Lagrange 
multiplier enforcing the same flatness condition, and the theory (for closed manifolds) amounts to an evaluation of 
the volume of the moduli space of flat connection, on the give topology Ai . 

The GFT framework (for the class of models we are going to consider) encodes the same quantum dynamics, by 
assigning amplitudes to any Feynman diagram dual to a given D-dimensional simplicial complex corresponding to 
the discrete path integral for the discrete version of the above action [17] ; on top of this, it embeds this dynamics 
in a wider framework which includes degrees of freedom coming from topology change (and singular complexes). No 
direct relation between the GFT equations of motion and the classical (continuum) BF equations of motion has been 
spelled out yet, though. 



A. Single-field GFT 

To begin with, let us consider what happens in the simplest case of a single field (non-colored GFT). 
The basic field is a function of a pair of group variables, 0(<7i,<72)) with g £ SU(2). The action we take is the one 
generating triangulated surfaces in the Feynman expansion for the partition function: 

S = - [ dg a dg b (/) ab 4> ba + - [ dg a dg b dg c (p ah 4> ca (l) bc , (42) 



where we use the notation cj) ab = <p(g a ,g b ) in order to avoid clutter. Note that we do not assume that cj) ab = 4>ba The 
equation of motion obtained upon variation of this action is: 



4>{h2, hJ + XJ dh 3 (f>(h 2 , ha)4>(h3,hi) = 0- (43) 

Notice also that if we had considered the interaction term containing, for example, the pairings 4> a b<pac4'bc, we would 
have obtained a different equation of motion. For a more detailed discussion of this point, see appendix [B] For the 
present discussion, we have chosen a pairing that ensures orientability of the simplicial complexes obtained in the 
GFT Feynman expansion. We consider now the heat kernel (coherent) states, in the 2D case, corresponding in the 
GFT context to bivalent vertices. Let us consider, as a trial configuration, 

t,G a ,G b {gi,9z]t) = i" J dhp t (gihG a )p t (g 2 hG b ), (44) 

which is obtained by the convolution of two heat kernels (with the same spread) peaked on different elements of the 
(complexified) Lie group, G^ 1 and G^ 1 (thus, in our case SL(2, C) elements). Here fi is a normalization parameter 
to be fixed by the equations of motion. Using the properties of the heat kernel, 

£(<?!, 52! t) = ^tigig^G^Ga). (45) 
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This latter is in fact a heat kernel on the coset SU(2) x SU(2)/L> 2 ~ SU(2V where D 2 = {(h,h)\h G SU(2)} C 
SU(2) x SU(2) (for a complete discussion of heat kernels on coset spaces, see [43|). 

It is clear that, if the two group elements G a ,Gb are the same, the field will be given by the heat kernel peaked 
around the identity in SL(2,C). Furthermore, it is obvious from this preliminary discussion that we will not be able 
to fix both G a and Gb- gauge invariance is telling us that one of them is redundant, i.e. the state will depend only 
on the equivalence class of elements under gauge transformations, labelled by G a G~^ 1 . However, we will be able to fix 
the value of the product G b ~ 1 G a using the GFT equations of motion. 

We now have to plug the trial solution (|44p into the equation of motion (|4"3"|) . For simplicity in the notation, we 
denote the product G^ G a = G. We get: 

A 4 [p2t{929^G) + X Wu (g 2 g^G 2 )] = 0. (46) 

Clearly, fj, — is a solution to this equation. However, this particular solution represents the "trivial" GFT vacuum 
associated to the absence of any geometry, and, actually, of any space at all; it is the trivial "Fock"GFT vacuum. 
We are interested, then, in the vanishing of the other term in the equation of motion. To analyze its content, define 

.9 = 9192 1 '- 

P2t {gG) + Xfip^gG 2 ) = V 5 e SU(2). (47) 

Being valid for any group element g, this equation corresponds actually to an infinite tower of equations (one 
for each value of g) for the classical phase space variables G, with parameters t, A and fi. These equations are 
then our "geometrodynamics" equations (in the simple BF case), obtained here as equations for the order parameters 
characterizing the GFT (mean field) background configuration in our (hydrodynamic) approximation, i.e. neglecting 
the contribution coming from fluctuations around the mean field, and coming directly from the full GFT dynamics. 



We are not able, at the present stage, to recast the above equations into a more geometrically transparent form, or 
to make more direct contact with the classical BF theory equations. However, we can try to solve them, and hope to 
make contact with classical BF theory at least at the level of the solutions we find. 

It is immediate to realize that an exact solution to this equation is t — 0, /i = — 1/A, G = I, i.e. 

6,i(ffi,ff 2 ;0) = -i% 1 ^ 1 ). (48) 

Despite its obvious relevance, this solution is not really admissible because it corresponds to an infinite total action 
(and thus lies outside the space of fields for which the GFT model is defined). Therefore, we keep t ^ and find 
some optimal value for the other free parameters of the heat kernel such that it is a reasonable approximation for a 
solution. In order to do so, we expand the above equation in representations, using the Peter- Weyl decomposition of 
the heat kernel: 

Pt(0) = '52d j e- tC 'Xj($), Xi(g)=tLn(9)S mn , Cj=j(j + 1). (49) 

3 

Using this decomposition in the equation of motion, we get a tower of equations (labeled by j) for t and G: 

tl m {G) + \pe- 2tc Hi m {G 2 )=Q, (50) 

or, equivalently, 

5 nm +Xpe- 2tc ni m {G) = {). (51) 

It is immediate to realize that the only SL(2,C) elements whose representation matrices are proportional to the 
identity for all representations j are those belonging to its center, namely I and —I. Indeed, their representation 
matrices are: 

ti m (±I) = (±l) 2j S nm . (52) 

The equation for j — provides the condition /i = —1/A, whence: 

1 - e-' 2tc i (±l) 2j = 0, j > 0. (53) 

Obviously, the exponential is always positive and smaller than one (since tCj > 0). Clearly, if G = —I, when 
specializing to the case of half integer j, the LHS of this equation would be always of order one, therefore failing to 
approach zero for any representation. Therefore, we have to keep G = I as the only plausible option. 



13 



However, unless t = 0, the GFT equations cannot be satisfied. In particular, for fixed i, for representations j is 
such that t\j > 1 we will have that the LHS of the equations are of order one. Therefore, our heat kernel with finite 
spread will represent an approximate solution to the equations of motion, as long as the spread is kept very small 
(mimicking, then, the exact solution given by the Dirac delta). The fact that this function is not a solution will be 
evident examining representations of sufficiently high spin, the threshold been given by the condition 2t\j 1. One 
has also to keep in mind that a spread of order t in the g (connection) variables implies a spread of order \jt in the 
conjugate B variables. 

With these considerations in mind, we can reach some partial conclusions. The trial field profiles (|44[) . with t « and 
G = I, can be seen as approximate solutions to the GFT field equation. Recall that G = I means G a = Gb in terms of 
the original GFT variables, which still implies that any wave function associated to a generic graph will have G = I as 
the phase space point associated to any graph link. They represent physical configurations of the (discrete) BF fields, 
that are sharply peaked around the trivial connection (everywhere on the graph), and peaked with large fluctuations 
a degenerate classical configuration of the B field, B — everywhere. Once more, we cannot make yet explicit 
contact between our effective "geometrodynamics" and classical BF theory, but one should notice that this selected 
configuration is one of the classical solutions of BF theory: a trivial connection, up to gauge transformations, and any 
B field satisfying B l (x)Bi(x) — const (the remaining components in the Lie algebra are pure gauge). So the candidate 
GFT configuration corresponding to LQG coherent states is a solution of the classical GFT dynamics if it corresponds 
to coherent states peaked on a degenerate (in the B), fiat (with respect to the connection oj) configuration 7 . Note 
also, once more, that, despite corresponding to a somewhat degenerate classical configuration, it does corresponds to 
a highly non-trivial quantum state, receiving contributions from arbitrary spin or group excitations, and in principle 
defined for arbitrarily complicated graphs, thus very far from the GFT free (Fock) vacuum. Still the result could be 
interpreted as suggesting that the full classical dynamics should be looked for in the dynamics of GFT perturbations 
around this degenerate man field configuration, rather than in the mean field equations themselves. 



B. Colored GFT 



The next case, in order of complexity, is represented by the colored version of the same model 8 , in which we have 
three fields, labelled by R, G, B, such that the action reads: 

S = \ J dg a dg b (^f n + tfafa + 0f b 0f Q ) + A J (dgftfafafa. (54) 



The equations of motion are: 



C + A / 4>bc4>cadgc = 0; <f>g + A / 4>t4>td9c = 0; (j>b a + A / ij&<j>°dg c = 0. (55) 



As in the previous case, it is immediate to check that: 

<t>ab = -ySigag^GR); 4>° b = ^j-S(g a g b ~ 1 G G ); <j)f lb = ^-5{g a g b l G B ) (56) 

is an exact solution, provided that <jro~go~b = —1 and three conditions for the parameters, namely Gb, = GgGb = I, 
and cyclic permutations. Using these conditions, it is easy to obtain 

Gr = ( r I, G g = CgI, G b = (bI (57) 

where C-RiCgjCb are signs satisfying the condition CrCgCb = 1- However, once more one cannot consider the config- 
urations as acceptable as such, because they lead to a divergent GFT action. 

As in the case of a noncolored model, we can then consider the regularized case of heat kernels with a finite spread: 

€b = tab = ^Mgag^Ga), (58) 



7 Recall the B field plays the role of a Lagrange multiplier in both continuum and discrete formulations, so that large fluctuations around 
its classical value do not imply dynamical instabilities. 

8 The addition of colors implies several constraints on the combinatorics of the resulting Feynman diagrams [4(1 [47J , and is important 
from the point of view of GFT symmetries 48] , but it corresponds to the same Feynman amplitudes and thus is expected to correspond 
to the same classical BF theory. 
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and similar expressions for the other colors. One gets: 

VRPt(gb9a lG n) + VgMb J Pt{gb9c 1 GG)pt{9c9a 1 GB)dg c = (59) 

After straightforward manipulations, this equation becomes: 

(j,RPt(gGR) + \[iGtJ.BP2t(gGBGG) = 0, (60) 

while the other equations can be obtained by even permutations of R, G, B labels. In terms of representations, we 
obtain the following three towers of equations: 

mSmn + XpafiBe-^ni^GBGaGx 1 ) = 0, (61) 
fi G S mn + XfJ-B PR.e tXj t{ nn (G rG bGq 1 ) = 0, (62) 

PsSmn + XpRPGe-^tUGaGRGB 1 ) = 0. (63) 

Once more, these are our effective "geometrodynamics" equations for the order parameters Gx labelling classical 
phase space points (of BF theory), obtained from the mean field approximation of the fundamental GFT dynamics. 
Again, taking first j = one gets: 

PR = — ; Pg = — ; Pb = — ; o R a G a B = -1 (64) 

Following the same reasoning of the scalar case, we get: 

GbGqGj^ = I GrGbGq 1 = I GqGrG = I, (65) 

which are the same that we have obtained in the case of the Dirac delta. Therefore, the structure of the solutions to 
the equations is the same, as far as the group elements associated to the peaks are concerned. 

The same comments on the trial configuration made in the case of the single-field model apply to this more 
complicated case. These heat kernels are not exact solutions, even though they can approach with arbitrary accuracy 
(t —> 0) an (unphysical) exact solution (but keep in mind the spread 1/t in the conjugate B variable). The geometrical 
interpretation of the (unique) solution found is similar, so we do not repeat it, with the only exception that we now 
might have that some of the parallel transports are peaked on the SL(2, C) matrix —I, which represents a full rotation 
of 2tt: in the case of half-integer representations, it leads to an overall multiplication by —1. 

IV. 3-DIMENSIONAL GFT 

We can extend the previous analysis to any higher dimensional case. In particular, here we focus on the 3- 
dimensional case, namely, for the non-colored case, the Boulatov GFT model for 3d BF (1st order gravity), in the 
Euclidean signature The Feynman amplitudes of this model are simplicial path integrals for the discrete 3d 

version of the action (|4"0)l . or equivalently correspond to the Ponzano-Regge spin foam model 0, 0, EH, [13], in terms 
of a discrete triad and a discrete gravity connection. At the canonical level, the phase space is parametrized by such 
discrete triad, an su(2) element, and the holonomy of the connection, an SU(2) element, for each edge of the graph 
by means of which any spin network state is defined (see [2, ; 3| as well as (50j). 

A. Single-field case 

To begin with, we consider the simple (non-colored) model with a single field. The field will depend on three group 
elements, projected using the diagonal action of the group 



</>(9i, 92, 9a) = dh<S>(gih,g 2 h,g 3 h), (66) 

JSU(2) 

in order to impose gauge invariance (closure of the triangle edge vectors in the metric representation). 
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The action is such that, in a perturbative expansion of a path integral, the Feynman diagrams are (dual to) oriented 
simplicial decompositions of three dimensional manifolds of arbitrary topology, and is given by: 

S = \ J (^flO 3 0123 0321 + J J {dgf '0123015604530426, (67) 

where the obvious notation has been employed. The reader is referred to appendix [B] for a careful discussion of the 
possible alternative definitions of the model, for what concerns the combinatorics of field arguments. Our particular 
choice (1671) . among those ensuring orientability, leads to the simplest equations of motion, reducing the number of 
terms to be manipulated. 

As in the 2d case, we will try to find (possibly, approximate) solutions in the form of the selected trial functions, 
dependent on a certain number of geometric parameters. The (highly nonlocal) field equation for the GFT will be 
then turned into equations for the geometric parameters, providing thus a sort of geometrodynamics. 

In our case, the trial function is the convolution of heat kernels: 



£,abc = V J dhp t (g a hG a )pt(gbhG b )pt(g c hG c ), (68) 

where G a ,G b ,G c are SL(2,C) elements encoding the geometrical properties of the mean configuration on which the 
state £ is peaked. In particular, in this three dimensional case, the relation between these group elements and the 
classical phase space of three dimensional geometry is transparent. In a polar decomposition of the SL(2, C) matrices, 
the SU(2) part will be associated to the mean value of the parallel transport around which the state fluctuates, while 
the positive Hermitian part will be associated to the average triads/fluxes. 

After simple manipulations, the equation of motion can be put into the form : 

dhpt (53 hG a ) pt (92 hG b ) p t ( g\ hG c ) + 



+ A/! 2 / (dh 3 )p t {gih 1 G a )pt{g2h2G b )pt(g 3 h.3G c )p2t{h 1 1 h 2 )p2t{h 2 1 h 3 )p 2t {h 3 1 h 1 ) 



= 



(69) 



Coming back to our concrete problem, as in the 2d case the choice /1 = gives an exact solution, associated to 
the "trivial" GFT vacuum. In order to find other solutions, we decompose the equation into a tower of equations by 
expanding it into representations. The result is: 
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(70) 

where the Wigner's 3j symbols have been introduced 10 . To simplify the form of the equation, we have defined the 
tensor 



(& 3 )^ ri (/ ll )^ 2r2 (^)^ r 3( /l 3)P2 t (/ l r 1 ^)p2 t (^^3)p2 t (/ l 3^l). 



(71) 



This tensor can be re-written as: 
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(72) 



It is worth mentioning that the expression in the interaction term would change dramatically if we were to change the interaction term 
in the original action, for instance by replacing 0123 by </>23l (an even permutation of the group variables, preserving the orientation 
properties). First of all, the equation of motion changes, involving more terms that have different combinatorial structures. Second, as 
a consequence of this first point, the various convolutions of heat kernels will be such that some (possibly all) of the heat kernels of the 
form P2t(h~ 1 hj) will be replaced by p2t(h~ hjG a ), where G a = G~, Gji. 

It is worth noticing the crucial way in which the particular choices for the various terms in the action affect the way in which the 
G a ,Gi,,Gc enter this equation. If we had chosen a kinetic term with the pairing 0i23<Ai23, we would make this equation totally 
insensitive to the SL(2, C) elements. In the case we are considering, the equation is insensitive to Gf, only. We do not have a clear 
understanding of this feature of the GFT dynamics. 
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where the 6j symbols have been introduced. 

One of the 3j symbols factorizes away, leaving us with 



Jl J2 ]3 

ri r 2 r 3 



) + V/O'i. h, h^t^GaGZ^rA^GcG- 1 ) ( jj J j ) = v./:../,../,. (73) 



Here, for convenience, we have defined: 



f(h,hJs;t) =Y,dud j5 d h e-^+^ { f I I V (74) 
1^ I H J5 36 J 



456 

These are the geometrodynamics equations for the order parameters derived from the GFT dynamics. Once more, 
we are not able yet to relate these equations directly with the classical BF equations, nor to give them a more 
geometrically explicit form. We then turn to the problem of identifying some solutions, which could instead be given 
a geometric interpretation. 

In order for the equation of motion to be satisfied, we need: 

ti nn {G a GZ 1 ) = 4 n {G a GZ 1 )5 mn Vj, (75) 

i.e. that the matrices are diagonal for any j. The properties of the representations we are playing with (see appendix 
ICj) lead then to the conclusion: 

G G C - X = h «£C\0 (76) 

Using this in the equation of motion, we can get rid of the second 3j symbol to get: 

1 + A/i 2 f(ji , 32 , 33 i *)<4 (G a G- 1 )c& 3 (G a G- 1 ) = OVji , h , h , m x , m 3 . (77) 

Therefore, we can argue that c^^GaGj 1 ) = c^, (GaG' 1 ), i.e. that the matrices are not only diagonal, but propor- 
tional to the identity matrix, in any representation. As we have already stated, the only SL(2,C) elements whose 
representation matrices are proportional to the identity in every representations are ±1. 
Consider ji, j2, J3 = 0. Then one gets: 

A^-p (78) 

which fixes the normalization of the field profile. 

From its expression, it is manifest that, if t £ M, / > 0. Therefore, we get 

1 - /(7i,j2,J 3 ;t)(±l) !yi (±l) % = 0, (79) 

that can be approximately be satisfied if the second term on the LHS is always negative, for all the representations. 
We get then that GaG' 1 = I, thus fixing the geometric part that can be controlled by the equation of motion (EOM). 

Therefore we obtain that the coherent state GFT configuration is a solution of the GFT dynamics for G a = G c 
and for arbitrary G , or equivalently, taking into account the gauge invariance properties of the field, we can conclude 
that the solutions are parametrized by a single SL(2, C) element G = GhG" 1 . The dynamics is therefore clearly much 
richer than in the 2d case already examined, as one would expect. The case G = I would correspond once more to a 
degenerate 3d geometry, but we now see that other solutions are allowed. One should now attempt to prove that the 
only freedom is in the phase space variables corresponding to the discrete triad, while the above configurations all 
correspond to flat connections, as one would expect from classical BF theory. We leave this for further investigations. 
Finally, conditions on t must be obtained by the analysis of the closeness of / to 1, similarly to what we have done 
in the case of 2D models. This requires a careful examination of the properties of the function /, that we also leave 
for future work. 



B. Colored model 



The investigations about the symmetries of GFTs [48[ have shown that, in order to implement within a GFT model 
a symmetry that corresponds to 3D simplicial diffeomorphisms, the colored generalization of the GFT model [4o| is 
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necessary, in which one has a multiplet of dynamical fields, each labeled by a color index. The quantum dynamics 
still generates random orientable 3D complexes of arbitrary topology, and it has been shown that the coloring also 
leads to the absence many singular configurations with respect to the usual models [H, IU [H[ ; this combinatorial 
properties of the resulting Feynman diagrams were in fact the original motivation for introducing coloring in the GFT 
framework. 

We consider the simplest incarnation of the theory, defined by the following action: 

S = \ J (dgf [0f 23 ^2l + ^23^21 + ^23^21 + 0123^32l] + A J (^^23^56^26 $43 ■ (80) 

The analysis of the equations of motion goes exactly as in the case of the single-field model. However, there are 
some crucial differences that make this model much more complicated, and the extraction of geometrical information 
significantly more involved. 

Concerning the trial functions, as in the 2d case, the SL(2, C) peaks are different color by color, at least in principle: 

Sabc =»x[ dhp t (g a hGX) Pt (g b hGX)p t (g c hG?), (81) 

where X — R,G, B, V. This implies that the structure of the tensor il that we have introduced in the previous section 
gets completely modified. Besides possessing an obvious color index, it depends in a non-trivial way on various 
products of SL(2,C) elements. For instance, the equation for the color R reads: 

Mii / dhp t (g 3 hG^)p t {g 2 hG§)p t {g a hGf) + Xp G p B p v [ (dhf [p t {g 1 h 1 G^) Pt {g 2 h 2 Gi)p t (g 3 h 3 Gl)x 



x p 2t {h^h 3 G v a {G® )- x )p n {h^h x G^{Gl)- x ) P%t {K{ x h % G B c (G^)- 1 )] = (82) 

In a decomposition in representations analogous to the one used to treat the scalar case, there are two obvious 
differences. First of all, the multiplication of the equations by the representation matrices of the elements (Gf*) -1 
does not eliminate any of the matrices contracted with the analogous of the tensor Q. Furthermore, the latter has a 
dependence on the SL(2,C) elements such that it cannot be easily rewritten in terms of 3j and 6j symbols. These 
complications prevent us, at the moment, from identifying the solutions of the equations. We can only stress the fact 
that the classical dynamics for the order parameters we have obtained seems richer than in the single-field case. 

However, it is easy to see that if we impose the constraint Gf = Gf , then the analysis of the solutions goes exactly 
in the same way of the single-field case, producing exactly the same results. 

The only difference involves the different normalizations. It is immediate to realize that: px = ^xp, with ax some 
signs, and p = |A| _1 ^ 2 . The signs are constrained to satisfy orogvbQv = — sign(A). 



V. EFFECTIVE DYNAMICS AROUND THE HEAT KERNEL CONFIGURATION 



In the previous sections we have considered some aspects of the mean field approximation applied to GFTs, at 
leading order, i.e. ignoring fluctuations around the mean field configurations. In order to understand its limits, 
and to start characterizing the full dynamics of the system around these non-perturbative vacua, it is important to 
elucidate the features of the effective theory for the small deviations away from the mean field. This is also needed to 
understand the stability of the configurations, and hence their viability as (nonperturbative) "ground states" for the 
GFT. Furthermore, in light of the possibility of expressing any GFT quantum dynamics in spin foam representation 
(which amounts to the perturbative expansion of the same in Feynman amplitudes written as functions of group 
representations), this analysis should clarify some of the features of the spin foam models that one would obtain 
around the new vacua and the differences with the "fundamental" ones. In the following, we will consider such 
effective dynamics in the same cases in which we have performed the analysis of the equations of motion and obtained 
conditions for the heat kernel GFT configurations. However, it should be clear, by now, that the analysis can be 
generalized to higher dimensional cases without any conceptual difficulty. 



A. 2D - single-field model 



The easiest case is the two dimensional simple GFT (without colors). We will use the splitting of the field into a 
mean field part (the heat kernel satisfying the mean field conditions p = — 1/A, G = I) and a fluctuation (p (satisfying 
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the same diagonal gauge invariance as the original field): 



<Pab = £,(g a ,9b) + <p{g a ,gb), £,{g a ,gb) = -jP2t(g a g b *)• 

The action for the field generates an effective action for tp ab which has the simple form: 



(83) 



- / dg a dg b (ip ab ) 2 + 



dg a dgb£,abfab + A / dg a dg b dg, 



tPabt;bc£,ca + fab^bc^ca + -^fab^bc^ca 



(84) 



In this effective action, we have a term, linear in ip ab , which would vanish if were a solution. However, as we have 
seen, ^ ab is not exactly a solution and hence it will give a (possibly small) effective contribution to the equations of 
motion. By taking care of the factors of A, it is easy to see that these terms are of order 1/A: in these sense they are 
nonpcrturbative. The value of t must be tuned in such a way that this term becomes negligible. We will consider this 
as done. 

The kinetic term is completed by a quadratic term in ip ab which is obtained by convolution of the interaction term 
of the original theory with a single copy of the background solution £ ab . This means, and this is a first crucial and 
generic feature of the effective theories we obtain, that this term brings the full microscopic dynamics of the underlying 
GFT model (as well as the properties of the particular background solution considered) into the effective dynamics 
of the perturbations, already at the level of the effective kinetic term 
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y a bVcd1C abcd {dg) 



(85) 



The latter has the explicit form: 



K 



abed 



S(g a g c l )${gb9 d *) - 2S(g b g c 1 )p2t(g a 1 gd) = 5{g a g c x )K9b9 d x ) 



dg e dg f V(g a , g b , g c , g d , g e ,g f ) p 2 t(g e X 9f) 



(86) 

Here, the disappearance of A is somehow fictitious. In fact, the conditions imposed on the parameters of the heat 
kernels (our effective "geometrodynamics" equations) are tying together G, t and A. Let us finally note that, in the case 
of two dimensional theories, the effective interaction term is just the same as the one present for the full theory, and 
hence its simplicial interpretation, even if not obviously useful or appropriate (given that now we may have already 
generated an effective spacetime by means of the choice of nonperturbative vacuum) is unchanged. 

It is useful to rewrite the action in terms of group representations, using the properties of the representation 
functions of SU(2) (see appendix [C| or, for a more complete treatment, [52]), and the consequent expression: 



ah — ^ ^ mA»(Mj )■ 



(87) 



The effective action reads: 

He 



Y J Y 3 e 

j mn nm j 

dj dj 



-2tC, 



yj yj i y3 yj v J 

mn nm ~ oj mn nr rm 

3dj 



Cj = j(j + 1) 



(88) 



We have an effective dynamics with a new, non-trivial propagator. In this simple case, this is the only formal 
modification with respect to the fundamental (microscopic) GFT action. Obviously, this will result in modified spin 
foam amplitudes. The nature of the microscopic GFT model as a tower of infinite unitary matrix models with 
dimension N = 2j + 1 is retained in the effective dynamics. 

The diagonal nature of the kinetic term allows to re-write the same action in terms of a trivial one, by performing 
the rescaling: 



yj _v yj — (l _ 2p- 2tC i) 1 / 2 Y j 
mn ' mn \ ^ / mi 



(89) 



The resulting action can be considered as coinciding with the microscopic one, but in terms of a new, renormalized 
(or better, momentum dependent) coupling constant: 



A^A cff (^(l-2e- 2 * c 0- 3/2 A 



(90) 
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The rescaled matrices can be seen as the components in representation space of a "quasiparticle" field: 

uj(g) = J2 z Ltim(9), (91) 

3 

related to the original fluctuation field by a nonlocal transformation: 

w(ff)= [T(g,h) V (h)dh, T(g,h) = '£d j (l-2e- 2tc ^ 2 x j (gh- 1 ). (92) 

3 

The amplitudes of the effective theory can be immediately derived by making use of the Feynman rules for the 
corresponding matrix model (see (HI). In particular, the amplitude associated with a 2D simplicial complex T, dual 
to a fat graph obtained from the spin j matrix model will be: 

A(T,j) = ±\* s djV = 1 x F (l - 2e-^)- 3F/2 df\ (93) 

where F is the number of faces (triangles) of the complex, x(r) is its Euler characteristic, and S is a symmetry factor. 

Stability The analysis of the scalar case, with its extreme simplicity, allows us to partially address the issue of the 
stability of the heat kernel solutions. In particular, it appears that there is an instability for low spins, given that 

1 - 2e _2tCj < 0, (94) 

showing that the coefficients in front of the kinetic term for such low spins is negative. Notice that this is even more 
the case the more one approaches the regime t ~ 0, which is the one in which £ is indeed a solution of the GFT 
equations. This shows that the heat kernel GFT configurations can be used as good approximations of exact solutions, 
but they are probably perturbatively unstable: the stability matrix has negative eigenvalues. If this is the case, it 
would be interesting to study in more detail the physics behind this instability, as it may signal the breakdown of the 
GFT hydrodynamic approximation, and indirectly give insights on the physical meaning of the same. 



B. 2D - colored model 



In the case of the colored model, the presence of several fields complicates the situation slightly. Indeed, the 
microscopic GFT interaction term induces a quadratic term characterized by a mixing between colors. 
The effective action written in group representations reads: 




2A 

+ 2e~ 2tAj (a B G J mn R^ r tl m (GB) + cx R B 3 mn G : ' nr t : ' rm {G R ) + a G Rm n Bl r t 3 rm (G G )) + — G 3 mn R 3 nr B J rm (95) 

It is a tower (labeled by the irreducible representations of SU(2)) of three-matrix models with a standard interaction 
term, accompanied by a number of terms (in which the background field £ enters explicitly) generating oscillations 
between the colors. These terms introduce a form of nontrivial dynamics even if we truncate the effective action 
to the lowest order in A. It is worth stressing that these effects are directly induced by the GFT interaction term 
due to the presence of a nontrivial background. In the simplicial, discrete spacetime interpretation of the resulting 
perturbative Feynman amplitudes (once more, not necessarily appropriate in this context) these terms is that they 
generate oscillations of Id simplicial spaces (or bivalent spin networks) in coloring, while keeping the topology fixed. 

The diagonalization of this oscillating kinetic term can, however, be performed, at least in the case G R = Gb = 
Gg = I (which is the case in which the heat kernel GFT configuration is a solution of the mean field dynamics), once 
it is recognized that the quadratic term can be rewritten as: 




(R 3 ,G 3 ,B 3 ) a B p 3 1 a R pj , (96) 
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where (3j = e~ 2tCj . 

The possibility, if not need, for such diagonalization in field space is the only new feature of the colored model with 
respect to the single-field GFT. 

One has to introduce new matrices, which are linear combinations of the representation matrices R 3 ,G J ,B J (the 
mapping being an orthogonal transformation, since the matrix in the kinetic term (|96|) is real and symmetric). In 
terms of these new matrices, the effective action will have a trivial kinetic term, at the price of complicating the cubic 
interaction term. The transformation needed to put the kinetic term in canonical form depends on the representations, 
through the coefficient /3j. Notice that for j — > oo, the transformation reduces to the identy. This is due to the fact that 
for large j the off-diagonal terms are exponentially suppressed. Therefore, the largest effect is on small representations, 
the crossover scale J being determined by the condition 2tCj « 1. 

The fact that the transformation of the matrices are j-dependent has important implications in determining the 
propagating modes of the effective model. Let us be schematic. We started from some fields: 

( f I ab=Yl X "^ t ^ 9a9 b 1 ^ X fnn= d 3 J V'abttnnigagb 1 ), (97) 

3 

and we have discovered that the normal (truly propagating) modes are not given in terms of X 1 but of 

Y a > j = M aI X IJ (98) 

* ran ±y± j ^*-mni W-V 

where M aI are the matrices diagonalizing the kinetic term (|96|) . The matrices Ys can be then associated to fields: 

u a {g) = F m»'4rM = / dh^(h)K aI {h,g), (99) 

j,n,m 

a nonlocal linear relation between the two kinds of fields, controlled by the kernel 

K aI (h,g) = £ djMflLnih^ig), (100) 

This is very similar to the relation between quasiparticles and atoms in the case of condensed matter systems. The 
effective quanta associated to the oscillations around the given (geometrical) GFT background are clearly collective 
modes of the fundamental GFT quanta. 

C. 3D model(s) 

The discussion of the effective action for the three dimensional case goes along the same lines of the two dimensional 
case. The effective action is given by: 



"Jabc 



. If S 2 S 

£ ° / 2\54>abc54>def 



^abcfdef 



1 / 5 3 S 



3! \54>abc54>def54>gh; 



(Pabcfdeffghi) + 7 (^PabcV } aef<fdbf ( Pdec) , (101) 
&=£ / 4 



where we are using a notation in which (...) denotes the integration with respect to the relevant group variables. 

We have already mentioned the term linear in ip, which is not exactly zero unless £ is a solution of the equation 
of motion. However, by tuning t we may control the effect of this term and make it negligible. Therefore the first 
nontrivial term in our expansion come from the second variation of the interaction term evaluated on £. 

In the end, the effective action of the model (neglecting the linear term) reads: 

Sett = o / ^123^321 + 77 / (VabcPaeftidbftidec + Vabc^dbf £,aef ^dec + ^Pabc^Pdec^aef^dbf) + (102) 



(103) 



2 / r — 2 



+ A J (dg)Vl23¥?156^42 6 £453 + ~ j {dgf ¥>123<£156<P426<^453- 
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The main points to notice are the following, all concerning the new non trivial kinetic term: 

1) due to the expansion of the theory around a nontrivial background configuration, the fundamental GFT in- 
teraction term percolates down to the quadratic and the cubic part of the effective action; also, each insertion of 
background field leads to a multiplication by |A| -1 / 2 , which means that the quadratic term is of order A ; 

2) the three new terms contributing to the kinetic term, in addition to the trivial one, are not equivalent: they 
correspond to three distinct combinatorial structures. Moreover, they corresponds to oscillations in representation 
labels along the lines of propagation. 

However, contrary to simpler 2D case, it is not yet clear whether it is possible to envisage a transformation, that 
would put the kinetic term in a canonical diagonal form. The nature of the oscillations, as the Feynman amplitudes 
resulting from the effective action, can be studied in representation space. One has to use the conditions on £: 

A < 0; £obc = pTjy^ J dhpt{g a hG a )p t (g b hG b )p t (g c hG a ), (104) 

with G ai Gb arbitrary elements of SL(2, C), and the expansion 

^123 = f n To 2 ) 4\nM)4lnM)4ln 3 (9^ (^5) 



11213 \ ri y 2 r g 
{j} ' 

Once more, the detailed analysis would not add much to the present discussion. 

We therefore restrict our attention only to the mixing terms in the quadratic part. Let us consider just one of them: 

A J (PabcfaefUb)idec(dg) 6 = A J (p a bcVaTJtdbf£dec{dg) 6 = (-l)x (106) 

r 3 r 5 r 6 ,n 2 n 3 n 5 n 6 • 

(107) 



where we have introduced 



W?l&ltn 2 n 3 n 5 n 6 = / {d 9 f {dhf (g,)^ (g e )t>ln 6 (<?/) X 

x Pt(gdhiG a )p t (g e h 1 G b )p t (g c h 1 G c )p t (g d h 2 G a )p t (g b h2G b )p t (g f h 2 G c )} , (108) 

to (slightly) simplify the notation. We might want to use G a — G c , and gauge invariance to fix G b = I. Orthonormality 
of the representations is expressed as: 
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dg a tl\ ni {9aWr\nM = ^-P^S rir J nirl4 (109) 



■Jl 



In turn, this implies that we can reduce the expression to: 



Eviii2j3 [ Jl J2 J3 \ yjij 5 je I Jl J5 J6 \ yirji;h333siti 
y ^nin 2 n 3 I r r r I y ^n 1 n 5 n 6 I I y ' r 2 r 3 r 5 r 6 ,n 2 n 3 n 5 n 6 i V iiv 7 



or 



EY313233 vh3&3e> ~u/h3 3 353<> /'111") 

where: 



UT3i;hj 3 j5j6 _ [ Jl J 2 J3 \ [ Jl J5 J6 \ w j 2 jaj5je (111) 

vv l:n 2 n 3 n s n 6 / t \ T\ T 2 ^3 / \ T\ T§ r§ I r 2r 3 rnr6,n 2 n 3 nsne- \- L - L ^J 

The structure for the other terms in the quadratic part of the effective action is the same, with differences only in the 
orderings/combinatorial structure. We could summarize the structure that we obtain into the expression: 



_ V313233 V343536 v-313233343536 (^^'^^ 
2 ys ~nin 2 n3^ v n4n 5 n 6 i^' nin2n3nin5n6 , \ ±±u / 
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FIG. 4: Basic propagator and vertex 
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FIG. 5: The three contributions to the kinetic term: the gray rectangles correspond to insertions of heat kernels. It is understood 
that a sum is performed over the lines that are beginning and ending at one heat kernel insertion. 



where 

iriihhHJBje —0 irj\h333*3sie — Olx,3l34X u/3i\3233js36 , xhisx M/-? 2 ">3x33 34. 36 j_X333eX T/f/ j3 iJihjih \ fii/) 
'^n 1 n 2 n 3 n i n 5 n e , ~ ' K 'n 1 n 2 n 3 n i n 5 n fi ^ \ u u nm4 vv l-n 2 n 3 n 5 n B ^ u u n 2 n 5 rr 2\n 1 n 3 n i n 6 u u n 3 n e w 3;rn.n 3 Tl4Ti5 J V ±±4 V 

and where 

Oirjihhjijsje — XjijeX A3235X X3334 X ( 1 1 c\ 

JK 'n 1 7i 2 n 3 n4n 5 n e u u nme, u u n 2 ns u u n 3 m \'~ 1 ~'- J ) 

is just the original kinetic term. 

The entire dependence on the structure of the interaction term and on the background solutions, i.e. t and the 
various group elements of SL(2, C), is encoded in the tensors W%, W 2 , W 3 . While the explicit structure of the quadratic 
term is not helping much, still the processes it encodes are rather clear. Indeed, we could write down the Feynman 
diagrams of the emergent theory at once, starting from the structure of the Feynman diagrams of the fundamental 
theory, constructed out of the basic propagator and vertex in figure 4. 

Apart from the symmetry factors, the various induced diagrams can be obtained very easily by heat kernels insertions 
on the external legs, in all the possible ways. Therefore, the three oscillatory contributions to the kinetic term will 
correspond to the diagrams in figure 

Nonetheless, let us stress once more that these peculiar contributions to the kinetic term are encoding a nontrivial 
dynamics, not associated to any change in the structure of the graph depicting the combinatorial structure of the 
pairings of field arguments, but only affecting its associated representation labels. 

This very same analysis applies to the case of colored 3D GFTs. Without going into the explicit calculations (very 
similar to the scalar case) , we can immediately say that one additional feature will have to be added. The percolation 
of the original interaction term (combining triangles with different colors into a tetrahedron) onto the effective kinetic 
term will result not only in an effective dynamics that at the lowest order (free theory) will change the representations 
labelling the lines of propagation, but also in oscillations in the colorings, in a way similar to the case of 2D colored 
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models. The quadratic part of the effective action reads: 

\ J ( d 9) 3 Wabcfcba + VabcVcba + VabcVcba + fabc^ba] + 

i \ ( [",„-R ,„G C B C V . ,R C G ,B C V . , R C G C B ,V , 
+ A / {dg) YPabcVaefZ,dbf£,dec + VabcZaefVdbfZdec + ( PabcLefUbf ( Pdec + 

iabc^aeffdbf^dec + £at>c i Paef€dbf i Pdec + ^abc^aefPdbfPdec] > 

where the £ x are determined as in the previous section. The percolation of the interaction term of the fundamental 
theory into the effective quadratic term, and hence in the effective propagator is even more evident. 

VI. FUNDAMENTAL VS EFFECTIVE SYMMETRIES 

In any field theory, the choice of the vacuum is crucial for the understanding of the symmetries of the effective 
theory describing the dynamics of the excitations around that vacuum. Viceversa, the behaviour under symmetry 
transformation gives important informations about the physical nature of the chosen vacuum itself. In particular, 
the ground state defines an order parameter (or a set of order parameters) determining the symmetry group of the 
truncated theory. Nonetheless, it is important to remember that the effective theory for fluctuations around the 
new ground state remains invariant under the original symmetry group: the difference is in the representation of the 
symmetry group, which becomes nonlinear. 

Let us clarify this point. Consider a field theory (e.g. a GFT) for some fields collectively denoted by $ (in- 
ternal/tensor/spinor indices are omitted). Assume that the theory is invariant under the action of a symmetry 
transformation: 

$ -> <f>' = [/$, S[U$] = S[<f>] (If 7) 

The Fock vacuum $ = is obviously invariant under the action of U. 

Imagine that the new (non-perturbative) ground state is now not invariant under the action of the symmetry, but, 
for example, only transforms covariantly: 

[/$ ^$o • (118) 

When we consider the effective action that we obtain expanding around $o 

$ = $0 + (119) 

we obtain that, under the action of the symmetry group: 

$' = U$ => <t>' = $' - $ = U<p + (U - I)$ • (120) 

Notice that here we have written exact expressions. We have never referred to the fact that <f> is small, in any sense, 
or that U is an infinitesimal transformation. Then, the effective action for cf> is not invariant under cp —> U<f>, but 
rather it is invariant under the transformation 

[7^+(Z7-I)$o • (121) 

This must be expected whenever an order parameter transforming under some representation of the symmetry 
group of the theory is introduced. 

In our GFT case, the heat kernels are parametrized by the group elements around which they are peaked. In turn, 
they determine, by the mechanism we have just outlined, the shape of the symmetries of the effective theory around 
the chosen vacuum. 

Let us discuss further this important point, in the case of 3d colored GFT, where the implementation of (simplicial) 
diffeomorphism symmetry has been recently unraveled [48| . It has been shown that the GFT action is left invariant by 
global (from the QFT point of view) transformations (forming a quantum group) of the GFT field, which correspond 
(at the level of the Feynman amplitudes, given by simplicial gravity path integrals, to 3d simplicial diffeomorphisms 
53]: in terms of simplicial complexes appearing in the perturbative expansion, they correspond to independent 
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translations of the vertices [48| . Their geometric meaning is manifest in the non-commutative metric representation 
ofGFTs[l7j|. 

These transformations, in the group representation, take the form: 

<Xsi,52,53) -> exp[Tr(ei25i5^ 1 ]exp[Tr(e235253 _1 )] ex P[ Tr ( e 3iff35r 1 )]^ l (5i 5 52,53) (122) 

where are su(2) Lie algebra elements, and the trace is taken in the fundamental representation (the phases 
are given by the non-commutative plane waves of [54| . coming in turn from the quantum group Fourier transform 
[EH [EH], used in [l7| to obtain the non-commutative metric representation of GFTs), and where we neglected the color 
labels, which, although crucial for the definition and understanding of the symmetry, are not relevant for our present 
exposition. 

As we have mentioned, the presence of a nontrivial vacuum would change the action of these transformations. In 
particular, we could consider the action onto the heat kernel configurations we have studied in this paper. 

Before studying the coherent GFT configuration, let us study the simpler case 

^(9U9293) = J dhS( 9l hG a )S{g 2 hG b )S{g 3 hG c ) . (123) 

The transformed field reads: 

^1 23 = exp( ? ;Tr[e 12 G- 1 G b ]) exp(iTr[e 23 G b - 1 G c ]) exp( i Tr[e 31 G- 1 G Q ])^ 12 3, (124) 
which is equal to "0123 if and only if: 

G a = Gb = G c (125) 

Therefore, the field configuration obtained by means of a convolution of Dirac deltas, in the general case, transforms 
covariantly: it gets multiplied by a phase. However, the phase turns out to be just the identity if all the the elementary 
distributions used to construct the field are peaked the same group elements. 

This is an instance of a more general fact. The group elements entering the field configuration (around which we 
expand the theory) as parameters are also influencing the properties under symmetry transformations, allowing, for 
instance, for different phases. We can have three different phases: 

I. G a = Gt = G c , and the field configuration is invariant under the full group of transformations (I122j) . This represent 

a diffeo-invariant solution. 

II. Two of the group elements are equal, e.g. G a = Gb- Therefore, the symmetry group is broken down to a smaller 

subgroup (generated by the only ei 2 ). We can interpret this case as a configuration possessing an isometry. 

III. G a 7^ Gb 7^ G c , which corresponds to the completely anisotropic configuration, and to a purely covariant field 
configuration. 

The above discussion, even if it refers to a field configuration which is unphysical, clarifies the logic that has to be 
followed in determining the symmetry properties of any particular GFT configuration one is interested in. In the case 
of heat kernels (which, we recall, become delta functions in the t = limit), it is immediate to see that, under such 
GFT diffeo transformations: 

1) the heat kernel (coherent state) configuration is not left invariant, and transforms in the general way (|122p in 
which any GFT field transforms; 

2) a a consequence, it is not transformed into another heat kernel peaked around the diffeo-transformed group 
element. This can be seen, for example, by taking the spread t to zero, and checking that the function obtained in 
this way is not a Dirac delta. 

Thus, on the one hand, as we expect from their interpretation as (second quantized) wave functions approximating 
a general continuum metric, the GFT coherent state configurations do not possess any special isometry, in the 
general case, and transform covariantly under GFT diffeos; the only case in which an invariance under diffeos is 
obtained corresponds to their t = limit, and to the special case G a = Gb = G c , which, as we had seen, can be 
interpreted as a solution of the classical GFT dynamics (somehow expected to encode also the projection onto diffeo 
invariant states, from the canonical point of view) corresponding, however, to degenerate geometries. On the other 
hand, the transformed configuration does not peak either on diffeo-transformed canonical phase space points, so it 
somehow does not seem reproduce faithfully the action of diffeomorphisms, one would expect from the LQG (canonical 
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quantum gravity point of view). For instance, this is at odds with the transformation properties of the canonical LQG 
"condensate" representation studied in [57l . l5q . 

Another interesting field configuration, similar to the above ones, that is worth discussing briefly is the exact solution 
of the GFT dynamics in 3d given by the field: 



identified first in [38j . It is easy to verify (see also [59] ) that this field configuration is invariant under a subset of GFT 
diffcos (|122l) (generated by 631) which, together with the additional rotation invariance that can also be identified at 
the GFT level [48|, [59| , form a deformed Poincare invariance corresponding to the Drinfeld double quantum group 
DSU(2). This confirms the interpretation of the same GFT configuration as a quantum flat space used in [38j for 
the interpretation of the effective field theory for (reduced) GFT perturbations as a matter field theory on a non- 
commutative fiat space (see also the analogue construction in the 4d Lorentzian case in 39], but confirms also the 
general scheme for classifying GFT phases in terms of their symmetry properties, that we outlined above. 

The sketchy analysis that we have reported here shows concretely how 

• the non-perturbative vacuum one could chose (on physical or mathematical grounds) as the mean field config- 
uration relevant for the hydrodynamic/continuum approximation of the GFT dynamics (and as the tentative 
description of the "quantum space condensate"), will, in general, not be invariant under diffeomorphisms; 

• in the resulting effective theory, the action of diffeomorphisms is not the linear one but rather is of the nonlinear 
type we described; 

• the effective theories can be classified in phases, by means of the residual symmetry, if any, of the configuration 
around which we expand the theory. 

These facts should then be taken into account when trying to derive an effective dynamics of geometry from the 
fundamental GFT dynamics, if we expect the effective theory to be given by (some modified form of) Einstein's 
General Rclatvity, which is characterized by its symmetry under diffeomorphisms. 

On top of this, we have to add another feature that we have uncovered in the discussion of the effective theories 
for GFT perturbations. The symmetry that we are considering here is expressed in terms of the field if which is not 
diagonalising the kinetic term of the effective field theory. The true effective symmetry should instead be identified 
from a combination of the effective symmetry group for the field ip and of the Bogoliubov-like transformation that 
diagonalizes the kinetic term (needed to get the true propagating modes), defining the physical field ip. 



In this paper we have made a first tentative step towards the extraction of an effective classical geometrodynamics 
from group field theory hydrodynamics, using mean field theory techniques as applied to quantum GFT. 

This meant first of all identifying a candidate macroscopic GFT configuration with characteristic order parameters 
endowed with a geometric interpretation. This was chosen as the vertex building block of loop quantum gravity 
semi-classical coherent states associated to arbitrary graphs |19| . The coherent states constructed out of this, and 
associated to graphs with N vertices are then interpreted, in the analogy with condensates, as N-particle states. 

We have then used this candidate GFT configuration as the mean field around which to expand in an hydrodynamic 
approximation of the microscopic GFT model, and obtained the relevant consistency hydrodynamic equations (the 
GFT analogue of the Gross-Pitaevski equations adapted to this reference field). These, in turn, become equations 
for the geometric order parameters (identifying points on the classical phase space of gravity/BF theory) and other 
constants entering the definition of the background GFT field. Although their relation with the classical spacetime 
theory behind the GFT models considered (2d and 3d BF theory) is not yet entirely clear, these equations can be 
interpreted as a form of classical geometrodynamics, here obtained directly from the GFT dynamics. 

The special nature of heat kernel coherent states is crucial. Therefore, the problem arising at this point would be 
to find a dynamical mechanism that selects these states among all the possible ones. This task requires considerations 
that are beyond the limited scope of this paper, involving the detailed examination of the full path integral (and a 
better understanding of GFT perturbative renormalization (37j). 



•0(51,52,53)= / dhSigih) f(g 2 h)6(g 3 h) = f(g2g 1 1 )S(g 3 g 1 1 ) 




(126) 



Summary, discussion and outlook 



20 



The only solutions to these equations that we have been able to identify seem to correspond to degenerate geome- 
tries/B variables and flat connections, with quantum uncertainty on the latter being very small and on the former 
being correspondingly very large. As we pointed out in the text, this result is compatible with the classical dynamics 
of BF theory. However, three points have to be made in this respect: 1) the complexity of the equations, and of the 
corresponding solutions (thus the classical geometries selected) seem to grow considerably with the dimension of the 
spacetime they refer to, as we would expect; 2) the addition of "coloring" seem to allow for further geometric content, 
that is however not easy to elucidate; 3) the full geometric content of a many-particle quantum state constructed 
out of our reference wavefunction, contrary to the case of ordinary condensates, involves also a gluing operation, 
depending on the graph one wants to reconstruct, on our reference wavefunctions associated to spin network vertices; 
this plays an important role in the computation of geometric operators, but it plays no role in our construction, thus 
making the geometric interpretation of the identified configurations trickier. 

It must also be stressed that even a background configuration peaked on degenerate geometries as the one we 
identified is a highly non-trivial state from the point of view of the fundamental theory. It corresponds to a possibly 
continuum space in which geometry is everywhere defined, even if degenerate, made out of a possibly infinite number 
of microscopic GFT quanta (spin network vertices), whose fundamental degrees of freedom are all excited (e.g. the 
states are obtained by a sum over arbitrary group representations). By way of contrast, the microscopic GFT vacuum 
is a no-space state composed of no GFT quanta at all, where all geometric operators are identically vanishing. 

Despite their approximate nature, these solutions still allows to gain some better understanding of the theory. 
Indeed, the classical action itself (and hence the equations derived from it) can be seen as a saddle point /stationary 
phase approximation to the full path integral that defines the quantum gravity model. In the language of quantum 
effective actions, the classical action is only the lowest order contribution to the formal semiclassical expansion, and 
therefore even exact solutions to the classical equation of motion will have the role of approximate solutions for the 
full theory. Consequently, it is more significant to find general conditions according to which the desired function 
approximates reasonably well a solution, rather than an exact solution itself. 

This is true even in the case of BEC. The mean field approximation (MFA) gives a good qualitative description of 
the phenomena associated to Bose-Einstein condensation, but in certain conditions, it simply fails, due to instabilities 
appearing in systems with large inhomogeneities or large time derivatives (42[. In other words, only for some kind of 
situations and for some kind of observables, a MFA is adequate in giving the correct semiquantitative results. 

In the case of GFT, we expect that, in the full quantum/statistical problem, the classical action will receive 
important corrections and hence that even exact solutions to the classical equation of motion will be able to capture 
only a portion of the physical features, and possibly only under certain conditions. In the next section, we will see 
explicit signals of the approximate nature of the theory in the appearance of some instabilities, when looking at the 
linearized dynamics around these semiclassical configurations. 

Another point to note is the following. It is not obvious, at this stage, that the entire classical geometrodynamics 
has to be looked for in the mean field equations for the background GFT field selected. It is well possible that some 
(if not all) continuum geometric degrees of freedom are encoded in the GFT fluctuations around a given background 
field, e.g. one corresponding to a continuum but degenerate geometry, and that the latter has to be carefully chosen 
(or constrained by the GFT dynamics) so to admit a geometric re-writing and understanding of the former. 

Last, we note that our mean field equations gave also constraints on the relative value of the semi-classical parameter 
t, entering the definition of the coherent states chosen as mean field, and the GFT coupling constant A. This coupling 
constant, although presently not well-understood, has been linked from various points of view to the cosmological 
constant [12H14I . |60( , and can be in general guessed to be related to the coupling constants of the emerging geometrical 
theory (be it General relativity of some modification of it). For example, in the case of 3d GFT, if this is something 
resembling 3d gravity, the obvious candidate would be the dimensionless quantity G^A. While this is obviously only a 
speculation at present, we should also note that the same type of relation has been obtained in analog gravity models 
based on Bose condensates, where the semiclassicality (small depletion factor) was related to small values of GjvA 
emerging in the Newtonian limit |35| . 

Beyond the hydrodynamic (mean field) equations for the geometric order parameters, we have also extracted the 
effective dynamics for GFT perturbations around the chosen background, and pointed out some of their general 
properties, in particular the way they encode the microscopic dynamics and the differences between its spin foam 
expression and the original spin foam model. Finally we have discussed at some length the role of symmetries, in 
particular the recently identified GFT diffeomorphism symmetry, in the emergent effective dynamics, as compared to 
that in the fundamental GFT model. 

Much more work remains to be done, both to develop and improve the results we have obtained, and to test if the 
path we are suggesting towards the solution of the problem of the continuum in quantum gravity is the correct one. 
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We mention here a few lines of further research. 

The first focus should be the identification of the relevant macroscopic ground state for GFT hydrodynamics. Our 
choice is the most natural one from the point of view of LQG and spin foam models, as it corresponds to the semi- 
classical states used in both contexts for approximating macroscopic classical geometries from the pre-geometric data 
labeling kinematical LQG and SF boundary states. However, already at this level other choices are certainly possible, 
and the above states can be criticized on the grounds of being purely kinematical and possibly unstable (and thus 
not truly semi-classical), when dynamics is taken into account. This seems to be confirmed by our analysis of their 
behaviour under the GFT dynamics. A different type of criticism comes from taking seriously the GFT framework 
as a second quantization of LQG (and thus of spin foam models) and the condensed matter analogy. If GFT are 
understood as a second quantized framework for spin network states or simplices [H|-[l5j], with LQG or simplicial 
wave functions in turn interpreted as many-particle states, then the semi-classical LQG states whose building block 
we used as background configuration for GFT dynamics in mean field expansion (and providing our geometric order 
parameters) are, in a sense, "too semiclassical" . They correspond, in fact, to quantum many-particle states such that 
each particle (spin network vertex or fundamental simplex) is individually semi-classical. Each vertex state is then 
given, as we have seen, by a product of wave functions associated to its links, and each of them is a semi-classical state 
for the corresponding quantum degrees of freedom. From the point of view of condensed matter, then, this means 
choosing as macroscopic configuration a state in which all the fundamental atoms are behaving semi-classically. In 
particular, this means that such states approximate very well a certain type of (kinematical) classical observables, 
but not those observable ("extensive" ones) that depend additively on the number of links and vertices of the spin 
network graph. Examples of such observables are most geometric observables like areas or volumes. It is well-possible, 
and even likely, that the origin of the correct classical dynamics of spacetime and the very emergence of continuum 
spacetime structures in Quantum Gravity is due to quantum properties of the underlying building blocks of quantum 
space, and thus captured by states for such building blocks that are collectively semi-classical, but individually highly 
quantum. Once more, an example of this type of behaviour is that of quantum liquids and of Bose condensates 
in particular. This example would suggest that the relevant vacuum state in the GFT context should be a second 
quantized coherent state, not a first quantized one (corresponding to LQG coherent states). This gives a further 
motivation for developing a second quantized picture and Fock structure for GFT states, which is the prerequisite for 
constructing such vacuum states. 

A detailed analysis is needed for extracting the full geometric of the GFT-induced equations for the order parameters 
(classical phase space data) we have obtained, and of the effective spin foam dynamics for GFT perturbations. As 
for the first, a first step would be to try to separate the equations for the SL(2, C) group elements into equations 
for their su(2) and SU(2) components, given that they have the classical interpretation of B field and connection 
for the underlying BF theories. Assuming that our ansatz for the relevant vacuum state is the correct one, the 
final goal would be to establish a complete dictionary the GFT dynamics in mean field approximation, i.e. the 
GFT hydrodynamics, and the emergent classical theory (including GR). We are quite far from this goal yet, but an 
important asset is probably going to be the (non-commutative) metric representation of GFT, introduced in [17|, that 
brings the geometric content of GFT dynamics to the forefront, both at the level of the action and of the Feynman 
amplitudes. In particular, for the purpose of extracting the geometric content of the GFT mean field equations, it will 
be useful to develop a (non-commutative) metric representation of the LQG coherent states and of our chosen GFT 
vacuum. Work on this is in progress. Similarly, the elucidation of the geometric content of the effective dynamics of 
GFT perturbations will benefit from such reformulation. Prior to this, however, it is important to study in more depth 
the physical interpretation of the perturbations themselves, to clarify if they carry gravitational/geometric degrees 
of freedom or if they should rather be interpreted as emergent matter |38l440| . Clearly, their interpretation depends 
strongly on the interpretation of the GFT vacuum chosen, which may induce, as in analog gravity condensed matter 
models [I3,[28|, Ho|, a background geometry to which the perturbations couple. 

Finally, it is clear that the analysis presented in this paper, and the above-mentioned further steps, should be 
carried out in the physically relevant case of four spacetime dimension, and thus for 4d gravity GFT models. This on 
the one hand requires a better understanding of them than the one we have at present, and on the other hand it will 
contribute to such better understanding and development. 

In any case, we believe that the first steps we have taken in this work indicate a path, from the microscopic 
picture of quantum space provided by group field theories (and other related approaches to Quantum Gravity) to the 
macroscopic description of the same consolidated in General Relativity, that is worth following further. 
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Appendix A: SL(2, C) matrix algebra 

For convenience, we recall some basic properties of the SL(2, C) matrices, and in particular the polar decomposition. 
Consider M £ SL(2, C). There is a unique decomposition of it in terms of a SU(2) matrix U and a positive hermitian 
matrix with unit determinant P, M — UP. Consider the matrix: 

Q = M^M (Al) 

It is a hermitian matrix, and hence it can be diagonalized via a unitary matrix W, such that: 

Q = WDW~\ (A2) 

with D a diagonal matrix. It is immediate to realize that this matrix is positive, i.e. that the eigenvalues are (strictly) 
positive. This is due to the fact that the (diagonal) matrix elements of Q can be related to norms of vectors once one 
recognizes that: 

(x\Q\x) = (x\M*M\x) = || M\x) || 2 . (A3) 

Define the square root of D to be the positive diagonal matrix C such that C 2 — D. With this matrix, define the 
square root P of the matrix Q to be the hermitian matrix 

P = WCW- 1 (A4) 

Given that M is nonsingular, also C must be nonsingular. Hence, it will admit an inverse, C _1 . Let us consider, 
then, the matrix: 

U = MP' 1 (A5) 

Let us prove that it is unitary. Consider: 

UU* = M(P" 1 ) 2 M t = M{M^M)- X M^ = I (A6) 

This proves that u <E U(2). By keeping track of the signs of the determinants, one finds that U <G SU(2). 

Being a Hermitian matrix with unit determinant, belonging to SX(2,C), the matrix P admits the following repre- 
sentation: 

P = cxp(eV l ), (A7) 

with <7j the Pauli matrices 11 and e* real numbers (Einstein convention is assumed). 
These matrices have a nice transformation property. Consider: 

a = exp(eV,), a' ~ waw^ , (A8) 

with w a SU(2) matrix. Then: 

a = w exp(eVj)ttr = exp(e l waiW^). (A9) 
However, since wuiW^ is a Hermitian, traceless matrix, it can be written as: 

waiW f = R J i (w)a :j , (A10) 
where R l j(w) are just the rotation coefficients in the Lie algebra. Hence: 

a' = exp(e'Vi), e H = R)(w)e j . (All) 



Notice that these are Hermitian matrices. Anti-Hermitian matrices are sometimes used in the LQG literature. 
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FIG. 6: Oriented tetrahedron, with labeling of the edges. 



Appendix B: The action for 3D GFT 



It is worth to spend some words on a technical point that has been mentioned only briefly, i. e. the construction of 
the interaction term needed for the construction of the GFT model for three dimensional theories. We are going to 
discuss the issue step by step, in order to clarify some ambiguities that arise from a detailed analysis of the theory. 

The starting point is an oriented tetrahedron. The four faces are oriented in such a way that their sides are ordered, 
respectively, 123, 156, 264 and 453 (see picture^]). The desired GFT must be adapted to this combinatorial structure: 
the field arguments must be ordered in the appropriate way to respect the induced orientation on the sides of each 
face, and the gluings of the faces to form the tetrahedron itself. 

In building the model, we need a field 0123, which is not invariant under generic permutations of the arguments. 
Indeed, if this were the case, it would be impossible to encode any sort of orientation of the faces, and the perturbative 
expansion would include a larger class of simplicial complexes than the orientable ones. For the moment, we will assume 
no invariance at all under any form of permutation of indices. We will comment later on the possibility of having 
invariance under even permutations {i.e. those respecting the orientation). 

The kinetic term, providing the gluing of the faces of adjacent tetrahedra, in the perturbative expansion, must be 
such that triangles are glued, with normals pointing in opposite directions. The natural candidate, then, is: 



- , 123V y 321, 



(Bl) 



where the second field has arguments obtained via an odd permutation. Even here, we could consider seemingly 
equivalent alternatives: 

0123 0321! 01230213! 01230132- 

It is worth mentioning what happens when varying the kinetic term, for each choice: 

' (rf.g) a 01230321 = 20 c6a (B2) 
(rf.g) 3 0123 0213 = 20 6ac (B3) 



'Jabc 



'I. be 



'I. be 



01230132 = 20 acb (B4) 



Therefore, the kinetic term alone will induce slightly different terms in the equation of motion, and hence a choice 
has to be done. We have chosen to use the paring 0123 0321- 

The interaction term suffers from the same kind of ambiguities. Indeed, we can construct several inequivalent 
interaction terms, corresponding, naively, to the same combinatorial structure of the tetrahedron, but with rather 
different properties. The first vertex we consider is: 



- , 123V y 156V ; 426V v 453 



{dgf. (B5) 



It is the one that we have considered in the paper. However, it is easy to see that we can generate a certain number 
of other operators by reshuffling the group elements in each field by means of an even permutation: 



n 23 <PtT 1 7T 5 7T 6 7r' 7T ' <?V 



W {dgf. (B6) 
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To see that these vertices lead, in general, to different systems, we examine the first variation of this term. 

S 



»abc 



{61 S 2 5^4>T tl . K5 „ e 07^ Tr^ 0«t4' + (B8) 
^i^^01230iriir 5 7r 6 07ri'7r£<' + ( B1 °) 

To uncover the structure of this term, we need to use the fact that the permutations form a group. We will need 
to use the inverse permutations, to be denoted by the letters w,w',w": 

uj{ir{a))=a (B12) 

Therefore: 

• The first contribution can be easily derived. It reads: 

4*Ha'^z'^e4^'n' A 'K' h , K' G 4 > 7r' / l' n '^' K c (B13) 

• K^tjne implies: 

5i = 9u( a y, 95 = g u (b)i 56 = 9uo{c), (B14) 
so that the second contribution reads 

0^(0)23^4^' '(w(c))07ri'ir" '(a»(6))irj' (B15) 

Comparing this two terms, it is immediate to realize that they are distinct, in general: g a , 9b, 9c can appear in different 
positions among the arguments of the three fields. To see this explicitly, consider the case in which tt', tt" are just the 
identity. We would get the two terms: 

<t>it a -K5TT$ 0466045c) 0w(a)23042cj(c)04cj(b)3- (B16) 

For this to be the same term, we would need that also tt is the identity. This is the case for the term we have used 
to analyze the heat kernels within GFT. 

This brief analysis shows how there is an ambiguity in the construction of the interaction term. However, this 
ambiguity is not only formal, but results in different equations of motion, i.e. in different physical content of the 
theory. Of course, this will influence the partition function obtained from different choices, and, most importantly, 
the correlation functions of the models, given that they will obey different Schwinger-Dyson equations. 

It is worth mentioning, however, that if we find a particular solution to the EOM with a particular choice of 
orderings, and this solution is invariant under general permutation of the group variables, then it will be a solution 
for the EOM obtained from any choice of orderings in the construction of the kinetic term and the interaction term. 



Appendix C: Useful facts about SU(2) representations 

It is convenient here to reproduce some basic facts about representations of SU(2). For more details we will refer to 
[52| . The representation matrices that we need, t l mn , are labeled by a non-negative half integer I, and the convention 
on the matrix indices is that they run from —I to +1. In terms of SL2c matrices they have the following expressions: 

= a -( m +")/^+ n y+ m ^2 C(l,m,n,j) \ , (CI) 



™1 7 8 J 
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where the complex entries of the matrices are constrained to obey aS — [3-f = 1, from the definition of SL2C, the 
coefficient appearing in the sum is defined to be 



^■^j)- .,y ( ' +M !w7r )l( 'i;!'!'"^ - < c2 » 

— m — ny.yt + n — + m — j)\ 

and, finally, the sum is over the integers between: 

M = max(0, rn + n), N — min(Z + m, I + n). (C3) 

These matrices reduce to the matrix representations of all the unitary, irreducible representations of SU (2) when we 
restrict the group elements to belong to this subgroup of SL2C. 

The definition allows us to give a simple derivation of a result that is often needed in the manipulations discussed 
in this paper. It is rather straightforward to derive the following results: 

tin ( ^ ~i ) = (-D l+m (-l) l+n t l mn ( ° £ ) , (C4) 

and 

4n ( 7 ^ ) = ( 7 8 ) ^ C5 ^ 

Combining these two results, we can prove that: 

Os) = (-l^-iy+V^Gr 1 ), (ce) 

for any element g € SLC. This is the key result that we need to prove: 

s£/(2) dff4n (9)Cn>(9) = ( m l m >) ( n ^) (C7) 

where we have defined the following object 
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